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ABSTRACT 

The  construction,  comparison  and  analysis  of  three  distinct  strain  gauge  balance 
calibration  matrix  models  with  various  orders  of  the  calibration  equations  was 
conducted.  The  aims  of  the  investigation  were  to  identify  the  accuracy  of  the  three 
different  calibration  matrix  models  and  to  analyse  their  behaviour  with  different  data 
optimisation  techniques.  A  computer  program  written  in  the  C  and  X/Motif 
programming  language  has  been  developed  to  analyse  the  matrix  models.  Two 
different  least  squares  metiiods  and  four  optimisation  techniques  have  been 
implemented  within  the  software.  The  accuracy  of  each  calibration  model  is  evaluated 
using  two  statistical  estimation  methods.  It  was  found  that  all  three  balance  calibration 
models  had  similar  behavior  in  terms  of  accuracy.  The  accuracy  of  the  equation  in 
estimating  the  loads  experienced  by  the  balance  increases  as  the  order  of  the  calibration 
equation  increases. 
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Executive  Summary 

Wind  tunnels  are  one  of  the  primary  sources  of  aerodynamic  data  for  aerospace 
research.  The  Australian  Defence  Science  and  Technology  Organisation  (DSTO) 
operates  two  major  wind  tunnels  at  the  Aeronautical  and  Maritime  Research 
Laboratory  (AMRL),  one  covers  the  low  speed  regime  and  the  other  covers  the 
transonic  speed  regime.  Results  obtained  from  wind  tunnel  tests  are  used  in  many 
areas,  such  as  aerodynamic  research,  aircraft  design,  and  validation  for  computational 
fluid  dynamics. 

Achieving  a  high  level  of  accuracy  in  wind  tunnel  test  results  is  essential.  Accuracy  of 
the  results  depends  on  many  factors,  such  as  the  data  acquisition  system  and  the  force 
and  moment  measurement  system.  At  AMRL,  the  primary  force  and  moment 
measurement  system  is  the  multi-component,  internally  moimted,  strain  gauge 
balance. 

A  strain  gauge  balance  must  be  calibrated  before  it  can  be  used  to  measure  forces  and 
moments  in  the  wind  tunnel.  The  aim  of  the  balance  calibration  is  to  obtain  a  set  of 
calibration  coefficients  which  enable  the  voltage  output  of  the  balance  to  be  converted 
into  the  corresponding  forces  and  moments.  There  are  many  ways  to  describe  the 
relationship  between  the  forces  and  moments,  and  voltage  output  for  a  particular 
balance.  Due  to  the  imperfection  of  balance  design  and  manufacturing,  and  the 
combined  loading  condition  during  wind  tunnel  testing,  second  order  and  above 
calibration  models  are  generally  used  to  accoimt  for  the  interaction  effect  between 
different  components  of  the  balance.  As  the  order  of  the  calibration  model  increases,  so 
too  does  the  complexity  of  the  matiiematical  expressions.  For  example,  a  general  third 
order  calibration  model  for  a  six  component  strain  gauge  balance  has  a  total  of  198 
calibration  coefficients. 

Three  different  balance  calibration  models  with  different  order  calibration  equations 
are  investigated  in  this  report.  In  addition,  various  calibration  data  optimisation 
techniques  are  applied  to  different  calibration  models.  A  computer  program  has  been 
written  to  provide  an  efficient  method  for  performing  the  comparison  and  analysis. 

One  of  the  main  findings  was  that,  of  the  15  balance  calibration  equations  used,  the 
order  84  coefficient  and  order  96  coefficient  equations  provide  a  more  accurate 
estimation  than  the  lower  order  calibration  equations  for  the  relationship  between 
voltage  output  and  applied  load. 
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Nomenclature 


Symbols  which  are  not  listed  here  are  defined  and  described  in  the  corresponding 
section  of  the  report. 


Symbol 

[C] 

""  A  “ 

c 

[H] 

Hi 

Hi,p 

[R] 

Ri 

X 

Y 

Z 

1 

m 

n 

SCi 

5i,p 

X 

seli 

11 


Description _ 

Calibration  coefficient  matrix 

Approximated  calibration  coefficient  matrix 

Applied  load  matrix 

Applied  load  reading  of  z*  component 

The  estimated  load/ moment  of  the  z*  component  for  the  p* 

calibration  data  point 

The  measured  load/ moment  of  the  z*  component  for  the  p* 

calibration  data  point 

Voltage  output  matrix 

Voltage  output  reading  of  z*  component 

Axial  force  component 

Side  force  component 

Normal  force  component 

Rolling  moment  component 

Pitching  moment  component 

Yawing  moment  component 

Standard  error  of  the  calculated  load /moment  (dimensionless) 

The  difference  between  the  estimated  and  measured  load /moment 

Chauvenet's  criterion 

Standard  error  (with  dimensional  unit) 

Mean 


Subscript 

i  =  l 

Axial  force  component 

i  =  2 

Side  force  component 

i  =  3 

Normal  force  component 

i  =  4 

Rolling  moment  component 

i  =  5 

Pitching  moment  component 

i  =  6 

Yawing  moment  component 

_ 

Index  for  the  calibration  data. 
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1  Introduction 

This  work  was  carried  out  as  part  of  the  wind  tunnel  infrastructure  program  at  the 
Aeronautical  and  Maritime  Research  Laboratory  (AMRL).  The  aims  of  the  investigation 
were  to  identify  the  accuracy  of  different  strain  gauge  balance  calibration  matrix 
models  and  the  characteristics  of  the  model  when  used  in  combination  with  various 
data  optimisation  techniques.  Based  on  the  results  from  the  investigation,  an  optimum 
calibration  model  is  recommended  for  use  in  both  the  low  speed  and  transonic  wind 
tunnel  facilities  at  AMRL.  The  software  written  for  this  investigation  will  be  integrated 
into  the  existing  data  acquisition  software  to  enable  real-time  conversion  between 
strain  gauge  balance  voltage  output  and  the  forces  and  moments  experienced  by  the 
balance  during  the  wind  tunnel  test. 

Investigations  into  three  distinct  strain  gauge  balance  calibration  models  with  different 
order  calibration  equations  were  conducted.  A  computer  program  written  in  the  C  and 
the  X/ Motif  computer  language  was  developed  for  the  analysis.  Two  least  squares 
methods,  four  data  optimisation  techniques,  and  two  statistical  estimations  have  been 
implemented  within  the  computer  program.  The  computer  program  generates  a 
calibration  matrix  by  calculating  the  calibration  coefficients.  Using  the  calculated 
calibration  matrix,  a  reverse  calibration  is  applied  by  the  program  to  obtain  the 
estimated  forces  and  moments.  The  accuracy  of  the  calibration  model  is  evaluated 
based  on  the  ability  of  the  calibration  matrix  to  estimate  the  forces  and  moments 
compared  with  the  measured  values.  The  standard  error  of  the  data  set  is  used  as  an 
indicator  of  accuracy  for  each  calibration  model  in  the  computer  program. 

All  three  calibration  models  display  very  similar  behaviour  in  terms  of  accuracy  for 
different  equation  orders.  As  the  order  of  the  calibration  equation  increases,  so  too 
does  the  accuracy  of  the  model.  This  is  because  the  higher  order  models  provide  a 
more  comprehensive  description  of  the  interaction  effect  between  the  balance's 
components.  Both  the  2"^  order  84  coefficient  calibration  equation  and  the  3’^'^  order  96 
coefficient  calibration  equation  achieved  a  significant  reduction  in  standard  error 
compared  with  titie  2"'^  order  27  coefficient  and  tiie  3'''^  order  33  coefficient  equations.  As 
the  equation  order  increases  to  the  fotirth  order,  the  additional  amoimt  of  interaction 
effect  aCcoimted  for  by  the  model  compared  with  3''*  order  equations  is  expected  to  be 
minimal.  Hence,  it  is  suggested  that  fourth  order  and  above  calibration  models  are  not 
necessary. 

The  results  indicate  that  high  interactions  between  balance  load  components  may  lead 
to  diverging  results  in  the  reverse  calibration  procedure.  This  is  because  the  strain 
gauge  balance  voltage  output  for  a  particular  loading  condition  may  represent  either 
positive  or  negative  loads  for  a  particular  component. 

In  general,  it  is  recommended  that  optimisation  techniques,  which  require  the 
elimination  of  calibration  data  points,  should  not  be  used.  This  is  because,  those  data 
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points  being  eliminated  actually  represent  the  physical  behaviour  of  the  balance  or  the 
data  acquisition  system. 

Due  to  the  non-linear  nature  of  calibration  data,  results  showed  that  the  calibration 
data,  which  covered  only  a  narrow  range  of  load,  led  to  an  increase  in  the  calibration 
model's  accuracy.  This  was  because  the  regression  model  was  more  effective  in 
modelling  a  narrower  band  of  data  with  a  higher  degree  of  linearity. 

Additionally,  this  work  showed  that  the  number  of  calibration  data  points  has  a 
significant  impact  on  the  values  of  the  calculated  calibration  coefficients.  If  insufficient 
calibration  data  is  provided,  the  least  squares  regression  methods  may  fail  to  obtain  a 
set  of  calibration  coefficients.  The  calibration  model  may  also  fail  to  produce  an 
accurate  estimation  of  the  measured  forces  and  moments,  or  in  some  cases,  it  may 
produce  diverging  results  in  the  reverse  calibration  procedure. 


2  Balance  Calibration  Models 


There  are  many  ways  in  which  a  balance  calibration  model  may  be  defined.  This  report 
concentrated  on  three  different  models,  each  model  being  distinct.  The  following 
balance  calibration  models  have  been  investigated: 

1.  [R]  =  [C][H] 

2.  [H]  =  [q[R] 

3.  [H]  =  [C][R-H]  (this  is  a  general  representation  of  this  calibration  model) 

In  order  to  compare  the  accuracy  of  these  models  extensively,  the  first,  second  and 
third  order  calibration  equations  of  these  models  were  investigated.  (A  complete  listing 
of  all  equations  for  the  three  balance  calibration  models  is  given  in  Appendix  A.) 

The  following  table  is  a  summary  of  the  three  balance  calibration  models  and  the 
corresponding  orders  of  the  calibration  equations  investigated  in  this  report. 


Balance  Calibration  Model 

Order  of 
Equation 

Number  of 
Components 

[R]  =  [C][H] 

[H]  =  [C][R] 

[H]  =  [C][R-H] 

1st 

5 

5  coefficients 

5  coefficients 

6 

6  coefficients 

6  coefficients 

6  coefficients 

2nd 

5 

20  coefficients 

20  coefficients 

liMMaiiiiiiM 

6 

27  coefficients 

27  coefficients 

27  coefficients 

84  coefficients 

84  coefficients 

84  coefficients 

6 

33  coefficients 

33  coefficients 

33  coefficients 

96  coefficients 

96  coefficients 

96  coefficients 

Table  1.  Summary  of  balance  calibration  models  and  orders  of  calibration  equations 
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2.1  MODEL  1:  [R]  =  [C][H] 

The  model  ctirrently  used  for  strain  gauge  balances  in  the  wind  tunnels  in  the  Air 
Operation  Division  (AOD)  at  AMRL  is  in  the  form  of: 

M=[cIh] 

This  equation  describes  the  physical  relationship  between  load  and  strain  gauge  output 
voltage,  ie.  the  strain  gauge  voltage  is  a  hmction  of  the  applied  load. 

2.1.1  First  order  equations 

Both  six  component  and  five  component  first  order  calibration  equations  are  modelled 
by  the  software.  At  AMRL,  six  component  strain  gauge  balances  are  primarily  used  to 
measure  the  aerodynamic  forces  and  moments  of  aircraft  and  missile  models,  and  five 
component  strain  gauge  balances  are  used  to  measure  aerodynamic  forces  and 
moments  of  stores  released  from  aircraft  in  the  transonic  wind  tunnel.  The  five 
component  strain  gauge  balance  does  not  measure  the  axial  force  component  (X). 

2.1.1.1  First  order,  six  component  equation:  6  coefficients 

The  first  order  equation  consists  of  six  terms  for  each  component,  and  each  equation 
corresponds  to  an  individual  component  of  the  balance. 


R,  =  C,,//,  +  C,,H,  +  C,3/f3  +  C,,H,  + 


where,  i  =  1,. .  .,6 


2.1.1.2  First  order,  five  component  equation:  5  coefficients 


The  five  component  balance  calibration  equation  consists  of  Y,  Z,  1,  m  and  n 
components. 


R.  —  C12H2  +  C-^H^  +  Ci^H^  + 


where,  i  =  2,...,6 


2.1.2  Second  order  equations 

Two  second  order  equations  were  investigated,  the  27  coefficient  equation  and  the  84 
coefficient  equation  for  a  6  component  balance.  Additionally,  a  20  coefficient  equation 
for  a  5  component  balance  was  also  investigated. 

2.1.2.1  Second  order,  six  component  equation:  21  coefficients 

The  second  order  equation  includes  the  addition  of  square  and  cross  product  terms  but 
it  does  not  include  the  cross  product  of  absolute  terms.  This  equation  has  a  total  of  27 
calibration  coefficients  for  each  component. 
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where,  i  = 


~  Q.l^l  ^1,2^2  . ^1,6^6 


+  +  . + 

+ . 


2.1.22  Second  order,  five  component  equation:  20  coefficients 

The  axial  force  component  is  not  considered  in  the  five  component  balance  calibration 
equation.  This  equation  includes  both  the  square  and  cross  product  terms  but  there  are 
no  absolute  cross  product  terms. 

~  ^i.2^2  ^1,3^3  . ^i,6^6 

22-^2  ■*‘^/,33'^3  . ■*■^1,66 -^6 

^/,23^2^3  Q,24^2^4  . ^i,S6^5^6 

where,  i  =  2,... ,6 


2.1.2.3  Second  order,  six  component  equation:  84  coefficients 

The  second  order  84  coefficient  equation  is  based  on  the  27  coefficient  second  order 
equation  and  includes  the  cross  product  of  absolute  terms. 

■*■^<•,2^2  . ■*■^1,6^6 


+  q|i|l^.|  +  Q|2lK|  + +C,|e||iye| 

+  +  . +  C,eeH/ 


H,H,\ 

’^^1. 2121-^2  ^2 

. ■'■^/,6|6|^6 

l^6| 

■^^(,12^1^2  ■*■^,13^1^3  . 

\H,H, 

. ■^^i,|56||^5 

+  <^U|2| 

H,\H, 

. 

He\ 

+  ^.Mli2 

Mu, 

. 

1^6 

where,  i  =  1,...,6 


2.1.3  Third  order  equations 

Two  different  third  order  equations  were  investigated,  the  33  coefficient  equation  and 
the  96  coefficient  equation. 
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2.13,2  Third  order,  six  component  equation:  96  coefficients 

The  6  component  third  order  96  coefficient  equation  consists  of  single,  square,  cubic, 
cross  product  and  absolute  terms. 


~  ■^Q,2^2  . ■*"Q,6^6 

+  C,,|,||«,|  +  C,j,|H,|  + . 


+  C,,„H,^+C,^H/+ . +  C,,«/// 


-^1 1^1 1  ^;,2|2|  -^2 1^2 1  • 

. 

■^^1,12^1^2  13^1^3  . 

+  C,s6HsH, 

^i,|l2|  l'^l-^2|  ^i,|l3|  |■^1^3K  • 

. +C,w|W,H.| 

^;,il2|  -^1 1-^2 1  “*■  ^;,i|3|  1-^3 1  "*■  • 

. 

1-^1 1^2 

. +  C,,M.I«5|H6 

111-^1  ■*'^i,222'^2  . ■*'^/,666'^6 

"*"^(,11111  -^1  ^(,|222|  -^2  . 

where,  i  =  1,...,6 
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2.2  MODEL  2:  [H]  =  [C][R1 

Instead  of  defining  [H]  as  the  independent  variable  in  the  calibration  model,  the 
following  mathematical  models  treat  [R]  (the  strain  gauge  balance  output  voltage)  as 
the  independent  variable  in  the  calibration  equation.  Therefore,  this  equation  implies 
that  the  load  is  a  function  of  the  strain  gauge  output  voltage. 

[h]=[cM 


2.2.1  First  order  equations 

22.1.1  First  order,  six  component  equation:  6  coefficients 


~  ^1,2^2  ^1,3^3  ■*"  Q,5^5  ^1,6^6 


where,  i=l,...,6 


2.2.1.2  First  order,  five  component  equation:  5  coefficients 

Hi  =  Ci  2^2  3^3  4-^4  5^5  6^6 

where,  i=2,. .  .,6 

2.2.2  Second  order  equations 

Two  second  order  type  equations  were  investigated,  the  27  coefficient  equation  and  the 
84  coefficient  equation  for  a  6  component  balance.  Additionally,  a  20  coefficient 
equation  for  a  5  component  balance  was  also  investigated. 

2.2.2.1  Second  order,  six  component  equation:  27  coefficients 

~  ^1,1^1  ■*■^,2^2  . ■*■^1,6^6 

^1,1 1^1  ^/,22^2  . ^1,66^6 

12^1^2  ■*■^1,13^1^3  . 56^5^6 

where,  i=l,...,6 

2.2.2.2  Second  order,  five  component  equation:  20  coefficients 

~  ^i,2^2  ^,3-^3  . ^1,6^6 

■*‘Q,22-^2  ■*■^1,33^3  . ■*'Q,66'^6 

^,23^2^3  ^1,24^2^4  . "*■  ^i,56^5^6 

where,  i=2,. .  .6 
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2.22.3  Second  order,  six  component  equation:  84  coefficients 

Hi  =  C, +  C,  2i?2  + . + 

+  1^1 1  +  ^i,|2|  1^2!  + . +  Ci,\6\  1^6 1 


+  +^i,22^2^  + . +  C';,66^6^ 


+  1^1 1  ^i,2|2|-^2|^2l  +  • 

. . +  <^;,6|6|'^6|^6 

■*“  ^U2^1^2  ^1,13^1^3  . 

■*■^(,56^5^6 

+  <^;,|12||^I^2|  +  C'/,|13||^1^3|  +  - 

. 

+  <^;,l|2|-^l|^2|  +  <^i,ll3|'^l|^3|  +  - 

. "^^/.slel^sl^el 

+  Q|1|2  1^1 1^2  +<^, -,1113  -^11^3  + 

. 1-^5 1-^6 

where,  i=l,...,6 

2.2.3  Third  order  equations 

Two  different  third  order  equations  were  investigated,  the  33  coefficient  equation  and 
the  96  coefficient  equation. 

2.2.3.1  Third  order,  six  component  equation:  33  coefficients 

Hi  —  +^(,2^2  ■*“ . ■*■^(,6^6 

■^^1,22^2  . "*“^1,66^6 

Q,I2^1^2  ^1,13^1^3  . 56^5^6 

■'■^i.222'^2  . 666-^6 

where,  i=l,...,6 
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1.132  Third  order,  six  component  equation:  96  coefficients 

~  ^1,1^1  ^1,2^2  . "^^,6^6 

1^1 1  +  ^,-,|2|  1^2 1  + . + C,M6|  1^6 1 


^i,22^2  . ^i,66^6 


+  ^/.1|1|  I^I  1  +  C'/.2|2|  ^2 1^2 1  +  ■ 

. 

^1,12^1^2  ^1,13^1^3  . 

.+  C,  56^5^6 

+  ^i,\\2\  1^1  ^^2 1  +  <^M13|  -^1^31  +  • 

. ■^^/,|56||'^5'^6| 

+  ^M|2i^l|^2|  +  C,.i|3|i?,|/?3|  +  . 

. +  <^,-,5|6|'^5|^6| 

+  ^i,|l|2KI^2  +Q|I|3|^1  ^3 

. '^^/,|5|6 1-^5 1-^6 

111-^1  222-^2  . ^1,666-^6 

+  ^/,|nil  -^1^  +C/.|222|  ^2^  + . +  ^/,|666|  ^6^ 

where,  i=l,...,6 

2.3  MODEL  3:  [H]  =  [C][R-H] 

This  calibration  model  assumes  the  loads  measured  by  the  balance  are  a  fimction  of 
both  balance  voltage  output  and  applied  load.  (lAI  Engineering  Division,  1998) 

Only  six  component  equations  are  considered  for  this  particular  balance  calibration 
model. 

2.3.1  First  order  equation 

~  ^i,2R^2  ^i,3R^3  ^i,5R^5  ^i,6R^6 

where,  i=l,...,6 

Note:  this  first  order  equation  is  the  same  as  the  6  component,  first  order  equation  of 
the  [H]=[C][R],  balance  calibration  model,  (see  Section  2.2.1.1) 

2.3.2  Second  order  equations 

Two  second  order  type  equations  were  investigated,  the  27  coefficient  equation  and  the 
84  coefficient  equation. 
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23.2.1  Second  order  equation:  27  coefficients 

Hi  —  ^i,2R^2  . "^^i.eR^e 

. +C,..H,'‘) 

+ . +  C,^HsHt) 

where,  i=l/-”/6 

2.3.2.2  Second  order  equation:  84  coefficients 

Hi  =  C„l/f^l  +<^i,2R^2  + . +  Ci,6R^6 

-(C,-.|,||H.|  +  C,j,||ff2|+ . +  C,j,||H.|) 

-(c,,„H,"+C,j3H,'+ . +  c,,«h/) 

. + 

-{C,„H,H,+C,„H,H,  + . + 

“  (^i.ji2|  1^1^21  . 

-  (C„|,H,|H,|  +  C,,|,|H,lH,l+ . +  C, 

-  + . + ) 

where,  i=l,...,6 

2.3.3  Third  order  equations 

Two  different  third  order  equations  were  investigated,  the  33  coefficient  equation  and 
the  96  coefficient  equation. 

2.3.3.1  Third  order  equation:  33  coefficients 

Hi  —  +  c,  2j}f?2  . Q,6/!^6 

-(c,nH,^  +  C,,^H^^+ . +C, 

-{C,nH,H,+C,j,H,H,+ . J 

-{C,mfl/+C,,^»2+ . +  C,,«6H/) 

where,  i=l,...,6 
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2.3.3.2  Third  order  equation:  96  coefficients 


where,  i=l,...,6 


H .  —  +C,  2x^2 


•■*■^1,6/? -^6 


-(c,-j,|I«.|+c,mI».|+ . +c,„|H«|) 

. +C,,„H/) 

-(C,,,|,|H,|H,|+C,,|^H,|H,|+ . +C„.|,|H.|/f,|) 

-  + . +  ) 


'Ml2|l 


|  +  Q|I3||^1^3|  +  - 

—  ■*■^.•,1561 1^5  ^6 1; 

!"*'^i,l|3|^l  1-^3 

+  C,|,|3|i?.|/f3+. 


+ . +C,,«H/) 


■  +  C,J5|6|«=K) 


•  +  ^,-,16661  ^ 


3  Calculation  of  Least  Squares  Calibration  Coefficients 

Calibration  coefficients,  [C],  are  a  set  of  constants,  which  are  used  to  describe  the 
loading  characteristics  of  a  strain  gauge  balance.  To  obtain  an  accurate  description  of 
the  balance,  an  adequate  number  of  data  points  are  required.  This  number  of  data 
points  is  largely  dependent  on  the  balance  calibration  equipment  available  to  an 
organisation.  The  distribution  of  applied  loads  should  cover  the  maximum  range  of  the 
balance  and  ideally  it  would  be  similar  to  the  loads  experienced  by  the  balance  during 
wind  tunnel  tests. 

Using  various  types  of  regression  models,  a  set  of  calibration  coefficients  may  be 
obtained  from  a  set  of  load  data.  In  this  report,  two  different  types  of  least  squares 
regression  methods  have  been  used  to  obtain  a  set  of  calibration  coefficients.  Least 
squares  regression  allows  all  six  components  of  the  balance  to  be  loaded 
simultaneously.  Hence,  the  interactions  among  various  components  are  accoimted  for 
in  the  set  of  calibration  coefficients.  Additionally,  the  balance  can  be  loaded  in  any 
particular  order,  hence  a  random  and  arguably  more  realistic  loading  matrix  can  be 
applied. 

The  two  regression  methods  are  described  below  for  the  3'^^  order  calibration  equation, 
[R]  =  [C][H].  The  same  methodology  applies  to  the  other  calibration  models  and 
equations. 
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3.1  Multivariable  Regression  Method 

The  mathematical  expression  of  the  3^^  order  equation,  [R]=[C][H]  definition  is  given  in 
Section  2.I.3.I. 

The  following  assumptions  have  been  made  in  this  regression  method: 

1.  random  error  is  assiuned  to  be  zero; 

2.  the  observed  values  of  the  independent  variable  (in  this  example,  the  value  of  [H]) 
are  measured  without  error.  AH  error  is  in  [R]. 

Since  the  balance  has  six  components  (X,  Y,  Z,  1,  m,  n),  six  expressions  are  used  to 
represent  each  component  of  the  balance.  The  entire  set  of  p  data  points  can  be 
expressed  using  matrix  notation,  where  p  is  the  number  of  data  points,  as: 


^2,1 

^3,1 

^4, 

^5.. 

^1,2 

^2,2 

-^3,2 

^4,2 

^5,2 

^6,2 

^2,3 

^3.3 

^4.3 

^5,3 

^6,3 

Ap 

^2. 

^S.P 

^2., 

^3.1 

-  H,; 

^1,2 

^3.2 

- 

[«]= 

^1,3 

^2.3 

^3.3 

•••  H,; 

1 

1 

- 

Note:  the  size  of  matrix  [H]  depends  on  the  complexity  of  the  calibration  model.  For 
example,  in  the  3'‘^  order,  6  component,  96  coefficient  equation,  [H]  will  be  a  (p  x  96) 
matrix. 

Each  component  of  the  balance  is  represented  by  33  'linear'  and  'non-linear'  calibration 
coefficients.  The  calibration  coefficients  calculated  using  this  least  squares  method  are 
only  an  approximation.  This  is  because  random  errors  are  expected  to  exist  among  the 
data  set  due  to  various  sources,  such  as  electro-magnetic  interference  (EMI),  random 
vibration  on  the  test  rig  during  the  calibration  process,  and  errors  induced  in  the  data 
acquisition  and  processing. 
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Cu 

c 

'-'2,1 

c 

'^3,1 

Q,. 

Q,! 

Q,i 

^1,2 

c 

'-'2,2 

c 

'-3,2 

^4,2 

Q,2 

Q,2 

c 

'^1,3 

c 

'-'2,3 

c 

'^3,3 

^4,3 

c 

'-5,3 

^6,3 

■^1,666 

^2,666 

Q,666 

^4,666 

c 

'^5,666 

Q,666 

The  multivariable  regression  has  the  following  expression  (Sprent,  1969): 

[c =([/?r[^]r[^rw 


3.2  Ramaswamy  Least  Squares  Method 

The  mathematical  expression  of  the  3’'^  order  equation,  [R]=[C][H]  is  given  in  Section 
2.I.3.I. 


The  Ramaswamy  method  (Lam,  1989)  states  that  the  calibration  coefficients  are  found 
when  the  residual  between  the  measured  strain  gauge  output  and  that  obtained  from 
the  calibration  equation  is  a  minimum.  This  can  be  expressed  as: 


+  • 


^1,666 


where  i  =  l,...,6. 


For  this  particular  3'''*  order  model,  there  are  33  coefficients  for  each  component  of  the 
balance  and  p  equations. 


+  2,p 

+  C12H2P  + 

+  p  +  Cj^H^p 


^i.2^2,p  ^i,3^3,p 


+  •' 

"'^^i,666^6.p 

‘'^^i,666^6,p 

\h2,p 

'■  +  ^i,666^6.p 

k.p 

+  -+C,«H6,/ -«,,]■  W.,/ 


=  0 
=  0 
=  0 

=  0 


By  putting  the  equations  above  into  matrix  notation,  the  balance  calibration  coefficient 
matrix,  [C]  can  then  be  calculated  as  follows: 


where. 
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Iff.,ff., 

Iff.,ff3,  - 

Iff.,ff6./ 

Iff.,ffa, 

Iffi,ff.,  - 

Iff3,ff., 

Iff3,ff., 

Iff.,ff«,' 

Iff.,’ff., 

Iff^/ff., 

Iff 

6,p  “6,^ 

T»i, 

^1. 

Iffi,ff., 

Iff.,*., 

Iff.,*.,  ■ 

Iff., 

Iff.,^., 

Iff.,*.,  - 

Iff.,*., 

Iff., 

Iff3,ff., 

Iff.,*., 

Iff.,*., 

.Iff«, 

Iff.,’ff., 

!«../*., 

Iff 

3  D 

6,p 

'Cu 

r  c 

'-'2,1  '-'3,1 

c  c 

'-'4,1  '^5,1 

Q,2 

c  c 

'-'2,2  '-'3,2 

^4,2  ^5,2 

^6,2 

[c]= 

Q.3 

r  c 

'^2,3  '-'3,3 

Q,3  Q,3 

^6,3 

_Q,666  Q,666  ^3,666  ^4,666  Q,666 

Q,666  _ 

3.3  Five  Component  Strain  Gauge  Balance  Calibration  Equations 

The  multivariable  and  Ramaswamy  least  squares  methods  are  also  applicable  to  the 
calibration  equations  for  a  five  component  balance.  However,  due  to  the  mathematical 
characteristics  of  both  of  these  least  squares  methods  which  require  matrix  inversion, 
the  axial  component  must  be  removed  from  both  the  applied  loads  matrix  [H]  and  the 
voltage  output  matrix  [R]  before  these  least  squares  methods  can  be  used. 


4  Balance  Reverse  Calibration 


The  balance  reverse  calibration  process  uses  the  derived  balance  calibration 
coefficients,  [C],  to  calculate  the  load  experienced  by  the  balance  based  on  the  strain 
gauge  voltage  output.  Different  balance  calibration  models  require  different  reverse 
calibration  procedures.  The  procedures  applied  to  the  three  calibration  models  in 
Section  2  are  given  in  the  following  four  sections. 
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4.1  Model  1:  [R]  =  [C][H] 

Due  to  the  nature  of  this  particular  model's  equation,  an  iterative  reverse  calibration 
method  is  required  (Galway,  1980;  Cook,  1959).  A  brief  summary  of  the  method 
outlined  by  Cook  (1959)  is  described  below. 

[f']=[ciIf]+[c2Ig] 

where,  [Cl]  is  the  linear  calibration  coefficient  matrix; 

[C2]  is  the  non-linear  calibration  coefficient  matrix; 

[F']  is  the  apparent  loads  matrix; 

[F]  is  the  true  loads  matrix; 

[G]  is  the  true  loads  pairs  matrix. 

Hence,  [f]=[Cl]''{lF']-[C2lG][ 

[F]=[Cir[F']+MG]  where,  [Z)]=-{cir  [C2]} 

In  the  first  iteration,  it  is  assumed  there  is  no  interaction  between  components  of  the 
balance,  so  that; 

[f]  =  [F,  ]=  [Cl]"*  [F']  “  Step  1 

In  the  second  iteration  and  onwards,  the  interactions  between  components  are  taken 
into  consideration  in  the  reverse  calibration  process.  The  true  loads  pairs  matrix  [Gi] 
can  be  calculated  using  the  [Fi]  matrix. 

[F]-[Fj=[F,]+[i)IC,]  -Step2 

For  further  iterations,  step  two  is  repeated, 

[F,]=[F,]+[dIgJ  -Step3 

In  general  form,  the  reverse  calibration  process  can  be  written  as, 

|[fJ=[f,]+[dIg..,] 


This  iterative  process  is  repeated  until  the  values  of  [F]  converge.  In  general,  a 
converged  solution  can  be  obtained  after  between  two  and  ten  iterations,  depending  on 
the  accuracy  specified  for  the  converged  values. 

4.2  Model  2:  [H]  =  [C][R] 

Unlike  the  other  two  calibration  models,  this  particular  model  does  not  require  an 
iterative  reverse  calibration  procedure.  Instead,  the  true  load  experienced  by  the 
balance  can  be  calculated  directly  by  multiplying  the  calibration  coefficient  matrix  [C] 
by  the  strain  gauge  voltage  output  matrix  [R]. 
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This  gives  certain  advantages  over  the  iterative  reverse  calibration  procedure.  In  tire 
case  of  an  iterative  calibration  procedure,  the  matrix  which  requires  inversion  will 
become  larger,  hence  reducing  the  efficiency  and  accuracy  of  the  results  due  to  the 
inherited  inaccuracy  within  the  matrix  inversion  routine.  If  tire  matrix  being  inverted  is 
a  singular  matrix,  ^e  reverse  iterative  calibration  procedure  carmot  be  achieved.  Most 
importantly,  a  non-iterative  reverse  calibration  procedure  eliminates  any  possibilities 
of  a  diverged  solution  (see  Section  7). 

4.3  Model3;[H]  =  [C][R-H] 

This  particular  calibration  model  requires  an  iterative  reverse  calibration  procedure 
because  the  true  load  [H]  is  a  function  of  both  voltage  and  true  load  [R-H].  The  reverse 
calibration  procedure  is  very  similar  to  the  method  described  in  section  4.1. 

A  simplified  version  of  this  particular  model  is  as  follow: 

[fl']=[ClI«]-[C2lJ/'] 

where,  [H]  is  the  true  load  matrix, 

[R]  is  the  voltage  output  matrix, 

[H”]  is  the  2"'^  order  or  true  load  pairs  matrix  (obtained  from  [H]), 

[Cl]  is  the  linear  calibration  coefficient  matrix, 

[C2]  is  the  non-linear  calibration  coefficient  matrix. 

In  the  first  iteration,  it  is  assumed  that  there  is  no  interaction  between  different 
components  of  the  balance,  hence; 

[H,]=[ClI«] 

In  the  second  iteration,  the  value  of  [H"]  is  obtained  from  [Hi]. 


In  general  form: 


[//.]=  [aM-[c2lH"..,] 


The  iterative  process  is  repeated  until  the  value  of  [H]  converges.  Depending  on  the 
level  of  accuracy,  in  general,  [H]  will  converge  after  between  two  and  ten  iterations. 


4.4  Five  Component  Balance  Calibration  Equations 

To  apply  the  reverse  calibration  method,  listed  in  Section  4.1  and  4.2,  to  the  five 
component  1®‘  and  2"^*  order  balance  calibration  equations  of  the  [R]=[C][H]  and 
[H]=[C][R]  models,  the  five  component  matrix,  used  in  the  regression  method,  must  be 
converted  into  a  six  component  matrix  by  adding  zeros  and  one  to  the  axial  force 
component. 
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The  following  example  demonstrates  a  six  component  calibration  coefficient  matrix  for 
a  five  component  balance  calibration  equation. 

Model:  [R]=[C][H],  order  five  component  calibration  coefficient  matrix. 


Cl 

C2 

C3 

C4 

C5 

C6 

Hv 

"l. 00000 

nooono 

.XljQQOQa.-.-. 

o.ooono 

_._.DJ10QQD_._. 

.._._llQQQQQ._.i 

X 

Hv 

0.00000 

6.7728k -3 

1.38834e-5 

1.03006e-3 

3.04363e-6 

-2.83162e-4 

Hz 

0.00000 

7.09958e-5 

3.24193e-3 

-5.75365e-4 

2.501 19e- 4 

9.79163e-6 

H, 

0.00000 

-3.18513e-5 

6.99459e-5 

9.30475e-2 

6.49542e-4 

2.84688^-4 

H. 

0.00000 

-4.44862e-6 

-1.64417e-5 

1.03809e-3 

3.97225e-2 

8.74806e-4 

H„ 

0.00000 

8.73207e-5 

-6.00392e-6 

-1.35077e-2 

1.18139^-4 

7.422104e-2 

I 


^  Five  component  calibration  coefficient  matrix  obtained 
from  least  squares  regression  method. 

Six  component  calibration  coefficient  matrix  is 
formed  by  adding  the  axial  force  component 


The  above  six  component  calibration  coefficient  matrix  can  now  be  used  directly  with 
the  reverse  calibration  methods  (see  Section  4.1  and  4.2)  to  calculate  the  estimated 
forces  and  moments  measmed  by  a  five  component  strain  gauge  balance. 

This  approach  for  transforming  a  five  component  to  a  six  component  calibration 
coefficients  matrix  may  be  applied  to  any  order  of  calibration  equation. 


5  Statistical  Analysis 

The  two  statistical  indicators  used  to  analyse  the  accuracy  of  each  calibration  model  are 
the  standard  error  and  the  coefficient  of  multiple  correlation. 

5.1  Standard  Error 

Standard  error  provides  an  indication  of  the  accuracy  (the  degree  of  dispersion)  of  the 
calculated  loads  and  moments  using  the  approximated  calibration  coefficient  matrix, 
[C],  as  compared  with  the  measured  values.  This  parameter,  given  in  Equation  1,  can  be 
used  as  a  benchmark  to  compare  the  accuracy  of  various  balance  calibration  models. 
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where. 


sel;  = 


Equation  1 


f  is  the  niimber  of  degrees  of  freedom  in  the  calibration  equations,  {the 
number  of  degrees  of  freedom  is  equal  to  the  number  of  calibration  co^cients 
per  component), 

N  is  the  total  nmnber  of  data  points  used  in  the  calibration  data  set, 

Ri,p  is  the  measured  p*  component  of  the  strain  gauge  balance  output, 

Rj  p  is  the  calculated  (approximated)  force  or  moment  component. 


seli  is  the  standard  error  with  dimension  [Newton  or  Newton.meter]. 


The  standard  error  calculated  using  the  above  formula  has  a  dimension  of  Newton  or 
Newton-meter.  To  convert  it  into  a  dimensionless  parameter,  the  standard  error 
calculated  from  Equation  1  is  divided  by  the  corresponding  balance  component's 
maximum  design  load.  Hence, 

_  Standard  Errorj  (sel; ) [N  or  N  •  m]  iqq^ 

Maximum  Design  Load  Range;  [n  or  N  •  m]  Equation  2 

Similar  to  the  definition  of  standard  deviation,  standard  error  can  be  used  to  describe 
the  distribution  of  data  points.  For  example,  Ise  represents  68%  of  the  measured  loads 
and  moments  values,  2se  represents  95%  of  the  total  measured  values  and  3se 
represents  99.7%  of  the  overall  measured  values. 


Figure  1.  Standard  error  distribution 

5.2  Coefficient  of  Multiple  Correlation 

The  coefficient  of  multiple  correlation,  r,-,  given  in  Equation  3,  indicates  how  well  the 
calibration  equation  describes  the  relationship  between  the  outputs  of  the  strain  gauge 
balance  and  each  component's  loads.  It  also  indicates  the  ability  of  the  calibration 
equation  to  estimate  the  load  measured  by  the  strain  gauge  balance. 
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Equation  3 


A 

where,  ^  is  the  estimated  value  of  load/ moment, 

H-  p  is  the  mean  of  the  component's  load/moment  for  all  the  calibration  data, 
H;  p  is  the  measured  value  of  load/ moment, 
n  is  the  coefficient  of  multiple  correlation 

The  range  of  n  is  from  0  to  1.  If  n  =  1,  it  means  the  correlation  between  tibe  calibration 
model  and  the  measured  calibration  data  is  perfect. 


6  Data  Optimisation 

The  aim  of  data  optimisation  is  to  improve  the  accuracy  (reduce  the  scatter)  of  the 
estimated  loads  and  moments  as  compared  with  the  measured  values.  The  results 
obtained  from  the  optimisation  process  should  have  a  high  level  of  practicality.  In 
other  words,  a  zero  or  near  zero  standard  error  can  be  a  meaningless  representation  of 
the  accuracy  of  the  calibration  model  and  near  zero  standard  errors  could  occur  if  only 
a  small  number  of  calibration  data  points  have  been  sampled. 

In  this  report,  four  different  optimisation  techniques  have  been  used  individually  or  in 
combination  to  investigate  the  overall  effect  of  the  results  of  various  calibration  models 
on  the  standard  errors. 

6.1  'Zero'  Data  Filter  Optimisation 

The  function  of  the  'zero'  data  filter  is  to  eliminate  any  values  close  to  zero  in  the 
calibration  data.  This  may  be  desirable  because  close  to  zero  data  points  may  be  due  to 
background  noise  instead  of  the  actual  loads  or  moments  applied  to  the  balance.  This 
filter  is  applied  to  the  calibration  data  before  the  calculation  of  the  calibration 
coefficients. 

6.2  Standard  Error  Optimisation 

The  standard  error  optimisation  identifies  potential  "outliers",  as  shown  in  Figure  2, 
from  the  calibration  data  based  on  the  standard  error  calculated  for  each  individual 
component  of  the  balance. 
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The  standard  error  for  each  component  of  the  balance  must  be  converted  from  a 

percentage  to  the  corresponding  dimensional  unit,  h. 

_  sCj  X  maximum  design  load; 

__ 


s,.,  = 


where. 


h  is  the  outlier  tolerance, 

5i,p  is  the  difference  between  the  estimated  and  measured  load  /  moment, 
Ht,p  is  the  estimated  load  /  moment, 

A 

Hj  p  is  the  measured  load  /  moment. 


The  condition  selected  for  an  outlier  is: 


In  the  computer  program  (see  Section  8),  if  an  outlier  is  foimd  in  a  line  of  data  (ie.  one 
calibration  point),  the  entire  data  line  is  eliminated  from  the  calibration  coefficient 


calculation.  The  definition  of  outliers  in 
Figure  2. 


Hi 

(Measured  values) 
-  ve  +  ve 


graphical  representation  is  shown  in 

•  Estimated  value  using  calibration 
equation  (fail  within  the  outlier 
tolerance) 

★  Estimated  value  using  calibration 
equation  (fall  outside  the  outlier 
tolerance) 


Range  of  outlier 
tolerance  {±bi) 


Figure  2.  Standard  error  optimisation  -  definition  of  outlier 

6.3  Chauvenet's  Criterion  Optimisation 

After  the  calibration  coefficients  are  calculated,  an  outlier  elimination  process  based  on 
Chauvenet's  Criterion  can  be  applied  to  the  calibration  data  set  to  reduce  the  standard 
error  of  the  results.  Chauvenet's  Criterion  detects  and  eliminates  potential  outliers 
calculated  by  the  calibration  coefficient  matrix  [C]  through  the  reverse  calibration 
process.  (AIAA,  1995) 
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The  Chauvenet's  Criterion  defines  potential  outliers  using  the  following  relations. 


= 


The  condition  selected  for  an  outlier  isj _ 

where,  T  is  the  Chauvenet's  Criterion,  which  can  be  calculated  by  the 

expression; 

T=XA[lnMJ 

1=0 

where,  Ao  =  0.720185  Ai  =  0.674947 

A2  =  -0.0771831  As  =  0.00733435 
A4  =  -0.00040635  As  =  0.00000916028 

The  above  expression  is  only  a  curve-fit  equation  for  r  using  Chauvenet's  Criterion  for 
N<833,333,  where  N  is  the  total  number  of  data  points  used  in  the  data  set. 

6.4  Optimised  Calibration  Matrix  With  Non-Optimised  Calibration 
Data 

In  this  optimisation  technique,  a  set  of  optimised  calibration  coefficients  are  obtained 
by  either  the  standard  error  optimisation  (see  Section  6.2)  or  the  Chauvenet's  Criterion 
optimisation  (see  Section  6.3).  A  reverse  calibration  process  is  then  carried  out  on  the 
optimised  calibration  matrix  with  the  original  calibration  data  set.  The  aim  of  this 
optimisation  technique  is  to  investigate  the  relationship  between  the  optimised 
calibration  matrix  and  the  entire  set  of  original  (non-optimised)  calibration  data.  In 
effect,  this  shows  how  well  the  optimised  calibration  matrix  represents  the  original  full 
data  set. 


7  Balance  Calibration  Models  Analysis 

The  different  balance  calibration  equations  have  been  investigated  using  the  computer 
program  -  CALIB  (refer  to  Section  8),  which  has  been  developed  by  the  author.  The  aim 
of  the  investigation  was  to  identify  the  accuracy  of  each  individual  model  based  on  its 
standard  errors  (see  Section  5.1). 

In  the  ideal  situation,  where  the  strain  gauge  balance  has  no  interaction  between 
different  components,  a  simple  first  order  balance  calibration  mathematical  equation 
can  accurately  convert  the  balance  voltage  output  to  its  corresponding  load. 
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In  reality,  due  to  the  imperfection  in  balance  design,  manufacture  and  deformation 
vmder  load,  a  certain  degree  of  interaction  between  components  of  the  balance  always 
exists.  Hence,  a  higher  order  balance  calibration  mathematical  model  is  reqmred  to 
account  for  the  component  interaction.  There  are  many  ways  in  which  a  balance 
calibration  model  may  be  defined.  This  report  concentrated  on  three  different  models, 
each  model  being  distinct  (see  Section  2). 

Two  sets  of  calibration  data  have  been  used  for  the  analysis.  One  set  consisted  of  1886 
data  points  obtained  for  the  six  component  Aerotech  strain  gauge  balance,  shown  in 
Figure  3,  now  being  used  at  AMRL.  Another  set  of  calibration  data  consisted  of  329 
data  points  for  the  Collins  six  component  strain  gauge  balance,  shown  in  Figure  4.  The 
aim  of  using  two  sets  of  data  of  significantly  different  size  is  to  investigate  the  effect  of 
the  size  of  die  calibration  data  on  a  particular  calibration  model  in  terms  of  accuracy. 
However,  care  must  be  exercised  in  drawing  definitive  conclusions  about  the  models 
from  only  two  data  sets.  A  complete  set  of  graphs  showing  the  standard  errors  for  both 
data  sets,  for  each  of  the  models,  is  provided  in  Appendix  D. 
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Figure  4.  Collins  six  component  strain  gauge  balance  used  in  the  AMRL  low  speed  wind  tunnel 

7.1  Number  of  Calibration  Data  Points 

An  adequate  number  of  calibration  data  points  are  required  in  order  to  calculate 
calibration  coefficients  using  least  squares  regression  methods.  For  example, 
computation  of  the  [H]=[C][R-H]  2"**  order  84  coefficient  equation  with  one  standard 
error  optimisation  failed  during  the  reverse  calibration  process  for  the  Aerotech 
balance  (1886  data  points),  because  some  of  the  estimated  values  calculated  using  the 
optimised  calibration  matrix  failed  to  converge  after  6  iterations.  This  behaviour  was 
caused  by  the  elimination  of  a  large  number  of  the  calibration  data  points  (83%  of  the 
total  number  of  data  points)  after  the  one  standard  error  optimisation  process  (see 
Section  6.2).  This  significant  reduction  in  the  nxtmber  of  calibration  data  points  caused 
the  least  squares  regression  to  fail  to  generate  an  accurate  description  of  the  balance 
behaviour.  However,  for  the  2se  and  3se  optimisation  case,  only  21.7%  and  5.2%  of  the 
data  have  been  rejected  and  the  solution  converges,  (see  Figure  12  in  Appendix  D.3) 

All  three  models  have  similar  standard  error  values  for  die  same  corresponding  order 
of  the  calibration  equation.  From  the  analysis,  there  is  no  firm  indication  that  one 
particular  model  is  superior  to  the  others  in  terms  of  accuracy  (see  Appendix  D).  With 
the  1886  data  set,  using  a  higher  order  equation,  such  as  the  order  84  coefficient 
equation  and  the  3'"'^  order  96  coefficient  equation,  the  [H]=[C][R-H]  and  [R]=[C][H] 
models  have  a  slightly  lower  standard  error  (=3.64x10"^%)  compared  with  the 
[H]=[C][R]  model. 

In  terms  of  the  reverse  calibration  process,  model  [H]=[C][R]  has  the  advantage  of 
using  a  simple  non-iterative  process,  which  eliminates  the  possibility  of  diverging 
reverse  calibration  results.  In  the  other  two  models,  if  the  interaction  between  different 
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components  of  the  balance  is  relatively  large,  there  is  a  possibility  that  an  iterative 
reverse  calibration  process  may  fail  to  converge.  This  may  occur  when  a  particular 
voltage  output  from  the  balance,  has  both  a  positive  and  a  negative  force  or  moment 
for  a  particular  component. 

7.2  Effect  of  Balance  Calibration  Equation  Order 

As  the  order  of  the  balance  calibration  equation  increases,  the  accuracy  of  the  model 
also  increases.  This  is  represented  by  a  reduction  in  standard  error  for  each  balance 
component.  For  example,  the  standard  errors  for  the  [R]=[C][H]  model  are  shown  in 
Table  2  and  Figure  5.  Since  each  component  of  the  balance  has  a  different  degree  of 
interaction,  the  standard  error  also  varies  between  components.  A  second  order 
definition  with  27  coefficients  has  a  significant  reduction  in  standard  error,  especially 
in  components  with  high  interaction  behaviour.  This  is  shown  by  a  reduction  of  0.385% 
(0.60054%  to  0.21562%)  in  the  standard  error  for  the  rolling  moment  component  of  the 
Aerotech  balance,  as  compared  with  the  first  order  6  coefficient  equation.  A  further 
0.075%  (0.21562%  to  0.14070%)  reduction  in  standard  error  for  the  rolling  moment 
component  is  achieved  by  using  the  second  order  84  coefficient  equation,  and  0.079% 
(0.21562%  to  0.13685%)  by  using  the  third  order  96  coefficient  equation. 


Standard  Error  [%] 

HX 

HY 

HZ 

HI 

Hm 

Hn 

6  coeff. 

0.43587 

0.16990 

0.13214 

0.60054 

0.07482 

0.14332 

0.25943 

27  coeff. 

0.08537 

0.21562 

0.06154 

0.09877 

0.10595 

Icjclj.ljjJ 

0.05654 

0.11171 

0.08477 

0.21083 

0.05641 

0.09307 

0.10222 

84  coeff. 

0.07138 

0.14070 

0.04497 

0.05999 

0.07490 

0.05069 

0.07810 

0.06969 

0.13685 

0.04379 

0.05900 

0.07302 

Table  2.  Standard  error  for  the  [R]=[C][H]  balance  calibration  model  with  1886  data  points 


The  standard  error  reduces  as  the  order  of  the  calibration  equation  increases.  As  shown 
in  Table  2,  there  is  only  an  average  of  0.002%  (0.07490%  to  0.07302%)  improvement  in 
standard  error  between  the  second  order  84  coefficient  and  the  third  order  96 
coefficient  equation.  Although  the  average  improvement  in  accuracy  is  low, 
components  with  a  relatively  high  degree  of  interaction,  such  as  the  rolling  moment 
component,  HI,  achieved  a  more  significant  reduction  in  standard  error  of  0.004% 
(0.14070%  to  0.13685%)  for  the  3'''^  order,  96  coefficient  equation  of  the  [R]=[C][H] 
model  compared  with  the  2"^  order,  84  coefficient  equation  of  the  same  model. 
Therefore,  4re  use  of  the  3''^  order  96  coefficient  calibration  equation  can  further 
improve  the  accuracy  in  estimating  the  load  experienced  by  the  balance,  in  particular 
for  components  with  a  high  degree  of  interaction. 

Similar  trends  for  the  other  models  are  given  in  Table  3  and  Table  4. 
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First  order  6  coefficient 
equation 

“■“Second  order  27  coefficient 
equation 

Second  order  84  coefficient 
equation 

“■“Third  order  33  coefficient 
equation 

“■“Third  order  96  coefficient 
equation 


Figure  5.  Balance  Calibration  Model:  [R]=[C][H]  with  1886  data  points 


Standard  Error  [%] 

HX 

HY 

HZ 

HI 

Hm 

Hn 

6  coeff. 

0.43583 

0.16990 

0.13213 

0.60045 

0.07482 

0.14332 

0.25941 

27  coeff. 

0.08537 

0.21537 

0.06149 

0.09848 

0.10599 

33  coeff. 

0.08473 

0.21040 

0.05645 

0.09321 

0.10227 

84  coeff. 

0.05187 

0.08055 

0.07129 

0.14107 

0.04508 

96  coeff. 

0.05079 

0.07794 

0.06952 

0.13771 

0.04389 

0.05986 

0.07329 

Table  3.  Standard  error  for  the  [H]=[C][R]  balance  calibration  model  with  1886  data  points 


Standard  Error  [%] 

HX 

HY 

HZ 

HI 

Hm 

Hn 

6  coeff. 

0.43583 

0.16990 

0.13213 

0.60045 

0.07482 

0.14332 

0.25941 

0.11579 

0.08537 

0.21562 

0.06154 

0.09877 

0.10595 

33  coeff. 

0.05653 

0.11171 

0.08477 

0.21083 

0.05641 

84  coeff. 

0.05182 

0.08045 

0.07142 

0.14086 

0.04499 

0.05990 

0.07491 

96  coeff. 

0.05058 

0.07812 

0.06966 

0.13702 

0.04385 

0.05893 

0.07303 

Table  4.  Standard  error  for  the  [H]=[C][R-H]  balance  calibration  model  with  1886  data  points 
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7.3  Balance  Calibration  Coefficients  Calculation 

Both  the  multivariable  regression  and  Ramaswamy  least  squares  method  produced  the 
same  calibration  coefficients.  Both  methods  require  matrix  inversion  at  some  stage  of 
the  process,  therefore,  these  methods  will  fail  if  the  matrix  being  inverted  is  a  singular 
matrix.  Out  of  the  two  regression  models,  the  multivariable  regression  is  easier  to 
program  than  the  Ramaswamy  least  squares  method. 

7.4  Data  Optimisation 

7.4.1  Standard  Error  Optimisation 

Using  the  1886  data  set,  a  reduction  in  standard  errors  is  achieved  by  applying  the 
standard  error  optimisation  process  (see  Section  6.2).  In  all  three  models,  a  significant 
reduction  in  standard  error  is  achieved  by  applying  Ise  optimisation.  The  increase  in 
the  model  accuracy  is  due  to  the  large  amoimt  of  data  being  excluded  from  the 
coefficient  calculation  process.  For  the  1886  data  set,  the  standard  error  optimisation 
process  leads  to  an  85%  rejection  of  data  from  the  original  calibration  data  set. 
Although  the  standard  errors  for  each  component  of  the  model  in  estimating  the 
reduced  data  set  are  reduced  significantly,  the  large  amoimt  of  data  being  removed 
from  the  original  data  set,  may  lead  to  a  model  that  does  not  actually  represent  the 
balance  behaviour.  This  means  that  the  accuracy  in  modelling  the  range  of  loads  may 
be  greatly  reduced,  and  the  calibration  matrix  may  not  be  a  'good  fit'  to  the  data. 

Although  both  2se  and  3se  optimisation  achieve  a  lower  standard  error  in  all 
components  with  less  data  points  being  excluded  from  the  original  data  set  (21.5%  for 
2se  and  4.4%  for  3se  rejection  of  data  from  the  original  data  set),  such  criteria,  in 
excluding  certain  data  points  from  the  original  data  set,  are  not  recommended  because 
those  points  being  eliminated  may  represent  the  actual  behaviour  of  the  balance. 

Wifii  the  329  calibration  point  data  set,  the  standard  error  optimisation  technique 
actually  increases  the  standard  error  of  some  components  of  the  balance  instead  of 
decreasing  it.  It  is  believed  that  this  is  because,  with  a  much  smaller  set  of  calibration 
data,  the  effect  of  removing  data  points  has  a  more  significant  effect  on  the  final  results 
compared  with  a  large  calibration  data  set,  such  as  the  1886  data  set. 

From  these  observations,  it  is  recommended  that  standard  error  optimisation  should 
not  be  used  to  increase  the  accuracy  of  a  particular  calibration  model  because  of  the 
elimination  of  data  points,  which  may  represent  the  actual  behaviour  of  the  balance. 

7.4.2  Linear  Segmentation  of  Balance  Load  Range 

The  accuracy  of  the  calibration  model  can  be  increased  significantly  if  it  is  applied 
separately  within  a  smaller  load  range.  For  example,  if  the  design  load  range  of  the 
drag  component  of  a  particular  balance  is  ±1000N,  instead  of  using  one  calibration 
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coefficient  matrix  to  describe  the  entire  load  range,  separate  individual  calibration 
matrices  could  be  used  for  reduced  load  ranges.  This  argument  is  supported  by  the 
analysis  using  the  329  data  set.  Standard  errors  obtained  from  the  329  data  set  are 
mudi  lower  &an  the  1886  data  set,  and  this  is  because  the  329  data  set  only  covers 
loadings  from  -240N  to  +240N  as  compared  with  the  1886  data  set,  which  covers  - 
lOOON  to  +1000N.  This  would  further  imply  that  the  size  of  each  smaller  load  range 
should  be  selected  based  on  the  best  degree  of  linearity  of  the  balance  loading 
characteristics  within  that  particular  load  range. 

As  shown  in  Figiure  6,  for  those  regions  where  the  balance's  loading  characteristic  is 
relatively  linear,  a  wider  load  range  can  be  selected.  In  the  case  of  non-linear  loading 
characteristics,  the  size  of  the  load  range  can  be  reduced  to  obtain  the  same  level  of 
accuracy  as  in  the  linear  region.  The  main  reason  to  support  this  linear  segmentation 
method,  is  because  least  squares  regression  methods  represent  a  set  of  data  points 
using  a  straight  line,  and  if  the  data  points  have  a  non-linear  characteristic  the  least 
squares  methods  are  not  able  to  describe  the  relationship  as  accurately. 


[C]o.20O  [C]200-500  [CJsOMSO  [C]800-900  [C]900-1000 

[C]65O«)0 


Figure  6.  Sub-dividing  the  design  load  range  of  a  balance  to  improve  calibration  model  accuracy 


7.4.3  'Zero'  Data  Filter  Optimisation 

By  eliminating  data  points  which  have  a  value  close  to  zero,  no  significant  effect  is 
observed  on  the  standard  errors  of  each  calibration  model.  In  fact,  the  practice  of 
eliminating  data  points  which  are  close  to  zero,  may  have  a  negative  effect  on  the 
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accuracy  of  the  calibration  model.  This  is  because  a  reading  close  to  zero  may  actually 
be  due  to  the  interaction  effect  of  the  balance's  components  and  not  backgroimd  noise. 

7.4.4  Chauvenet's  Criterion  Optimisation 

In  comparison  with  the  standard  error  optimisation  procedure  (see  Section  6.2  and 
Section  7.4.1),  Chauvenet's  Criterion  eliminates  an  average  of  2.2%  of  data  points  from 
each  original  calibration  data  set  while  achieving  a  significant  reduction  in  the 
standard  errors  for  each  component  of  the  balance.  The  maximum  standard  error 
reduction  of  20.5%  is  achieved  in  the  order  6  coefficient  equation  as  shown  in 
Figure  7.  (Refer  to  Appendix  D  for  the  results  of  the  Chauvenet's  Criterion  for  various 
order  calibration  equations.) 

0.7 

0.6 

0.5 

i 

I  0.4 
III 

■s 

•§  0.3 
0.2 

0.1 

0 

HX  HY  HZ  HI  Hm  Hn 
Balance's  Component 

Figure  7.  Effect  of  Chauvenet's  Criterion  on  Standard  Errors  for  [R]=[C][H],  1*'  order  6 
coefficients  equation  with  1886  data  points 

Data  points  eliminated  by  Chauvenet's  Criterion  are  purely  based  on  statistical 
analysis  with  no  consideration  of  the  data  point's  representation  of  the  actual 
behaviour  of  the  balance.  Hence,  this  optimisation  technique  should  be  used  with  care. 
It  is  recommended  that  data  points  removed  by  Chauvenet's  Criterion  should  be 
documented  and  reviewed  manually  to  check  if  some  kind  of  physical  and/or 
theoretical  correlation  is  evident. 
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7.4.5  Optimised  calibration  matrix  with  non-optimised  calibration  data 

Using  a  standard  error  optimised  calibration  matrix  with  non-optimised  calibration 
data  for  the  reverse  calibration  process  resulted  in  a  lower  accuracy  (increase  in 
standard  errors)  compared  with  the  standard  error  optimised  calibration  matrix  with 
optimised  calibration  data.  This  behaviour  is  due  to  the  optimised  calibration  matrix 
only  being  able  to  accurately  estimate  the  loads  and  moments  for  the  optimised 
calibration  data.  If  calibration  data,  other  than  those  within  the  range  of  the  optimised 
data,  is  used,  the  optimised  calibration  matrix  produces  inaccurate  load  estimations. 

Figiue  8  and  Table  5  show  the  increase  in  standard  errors  when  loads  and  moments  are 
calculated  from  the  optimised  (Ise  optimisation)  calibration  matrix  with  non-optimised 
calibration  data  for  the  [H]=[C][R],  2"^  order  27  coefficient  equation.  Similar  behavioiu 
is  also  found  in  the  other  calibration  models  and  equations,  as  shown  in  Table  6  and 
Table  7.  Standard  errors  from  this  optimisation  technique  increase  significantly  for  all 
balance  components. 


^Optimised  calibration  matrix 
with  non-optimised  calibration 
data 


■  Non-optimised  calibration 
matrix  with  non-optimised 
caiibration  data 


■  Optimised  calibration  matrix 
with  optimised  calibration  data 


Balance  Component 

Figure  8.  Calibration  model:  [H]=[C][R],  2”^  order  ecjuations  with  1886  data  points 
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Calibration  Matrix 

Ise  Optimised  Calibration 
Data 

Non-optimised  Calibration 
Data 

Ise  optimised 
calibration  matrix 

sei 

0.03206% 

sei 

0.06598% 

se2 

0.05795% 

se2 

0.16158% 

ses 

0.04057% 

ses 

0.10531% 

se4 

0.11129% 

se4 

0.35291% 

ses 

0.02717% 

ses 

0.15965% 

see 

0.04546% 

see 

0.10913% 

Non-optimised 
calibration  matrix 

sei 

0.05858% 

se2 

0.11579% 

ses 

0.08537% 

se4 

0.21562% 

ses 

0.06154% 

see 

0.09877% 

Table  5.  Standard  errors  for  Ise  optimised  calibration  matrix  with  non-optimised  calibration 
data  for  the  [R]=[C][H],  2”‘^  order  calibration  equation  with  1886  data  points 


Calibration  Matrix 

Ise  Optimised  Calibration 
Data 

Non-optimised  Calibration 
Data 

Ise  optimised 
calibration  matrix 

sei 

0.03264% 

sei 

0.06421% 

se2 

0.05706% 

se2 

0.15216% 

ses 

0.04034% 

ses 

0.10593% 

se4 

0.11079% 

se4 

0.34769% 

ses 

0.02618% 

ses 

0.16481% 

see 

0.04572% 

see 

0.11289% 

Non-optimised 
calibration  matrix 

Bi 

sei 

0.05860% 

se2 

0.11663% 

ses 

0.08537% 

se4 

0.21537% 

ses 

0.06149% 

see 

0.09848% 

Table  6.  Standard  errors  for  Ise  optimised  calibration  matrix  with  non-optimised  calibration 
data  for  the  [H]=[C][R],  2’^  order  calibration  equation  with  1886  data  points 
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Calibration  Matrix 

Ise  Optimised  Calibration 
Data 

Non-optimised  Calibration 
Data 

Ise  optimised 
calibration  matrix 

sei 

0.03247% 

sei 

0.06639% 

se2 

0.05742% 

se2 

0.15661% 

sea 

0.04042% 

sea 

0.10463% 

se4 

0.11101% 

se4 

0.34818% 

ses 

0.02667% 

ses 

0.16417% 

see 

0.04556% 

see 

0.11018% 

Non-optimised 
calibration  matrix 

sei 

0.05858% 

se2 

0.11579% 

sea 

0.08537% 

se4 

0.21562% 

ses 

0.06154% 

see 

0.09877% 

Table  7.  Standard  errors  for  Ise  optimised  calibration  matrix  with  non-optimised  calibration 
data  for  the  [H]=[C][R-H],  order  calibration  equation  with  1886  data  points 


For  the  Chauvenef  s  Criterion  optimisation  technique,  no  significant  improvement  in 
accuracy  was  foxmd  for  the  optimised  calibration  matrix  with  non-optimised 
calibration  data  (see  Figure  9),  compared  with  the  non-optimised  calibration  results. 
This  is  because  the  Chauvenet's  Criterion  optimised  calibration  matrix  estimates  the 
loads  only  within  the  optimised  calibration  data  set,  and  because  very  few  points  are 
eliminated,  the  optimised  load  data  set  is  similar  to  the  full  data  set.  This  behaviour  can 
also  be  observed  for  other  calibration  models,  as  shown  in  Table  8. 


2 
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0.25 

0.2 

0.15 
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0.05 
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^Optimised  calibration  matrix 
with  non-optimised  calibration 
data 

® Non-optimised  caiibration 
matrix  with  non-optimised 
caiibration  data 

■  Optimised  calibration  matrix 
with  optimised  calibration  data 


HX  HY  HZ  Hi  Hm  Hn 


Balance  Component 


Figures.  Chauvenet's  Criterion  optimisation  for  calibration  model:  [H]=[C][R],  2’^  order 
equations  with  1886  data  points 
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Calibration  Matrix 

Chauvenet's  Criterion 
Optimised  Calibration  Data 

Non-optimised  Calibration 
Data 

IR]= 

:1CI[H1,  Z"**  order,  27  coefficient 

equation  | 

Chauvenet's  Criterion 
optimised  calibration 
matrix 

sei 

0.05862% 

sei 

0.05858% 

se2 

0.11711% 

se2 

0.11579% 

sea 

0.08546% 

ses 

0.08537% 

se4 

0.21672% 

se4 

0.21562% 

ses 

0.06168% 

ses 

0.06154% 

see 

0.09957% 

see 

0.09877% 

[HI: 

=[C]|R],  Z”'*  order,  27  coefficient  equation 

Chauvenet's  Criterion 
optimised  calibration 
matrix 

sei 

0.05866% 

sei 

0.05860% 

se2 

0.11798% 

se2 

0.11663% 

ses 

0.08548% 

ses 

0.08537% 

se4 

0.21654% 

se4 

0.21537% 

ses 

0.06161% 

ses 

0.06149% 

see 

0.09918% 

see 

0.09848% 

[H]=| 

CllR-H],  2"*'  order,  27  coefficient  equation 

Chauvenet's  Criterion 
optimised  calibration 
matrix 

sei 

0.05862% 

sei 

0.05858% 

se2 

0.11713% 

se2 

0.11579% 

ses 

0.08546% 

ses 

0.08537% 

se4 

0.21674% 

se4 

0.21562% 

ses 

0.06168% 

ses 

0.06154% 

see 

0.09956% 

see 

0.09877% 

Table  8.  Standard  errors  for  Chauvenet's  Criterion  optimised  calibration  matrix  with  optimised 
and  non-optimised  calibration  data  for  the  2"^  order  calibration  equation  with  1886  data 
points 


In  summary,  it  is  recommended  that  the  optimisation  techniques  described  in  this 
report  should  be  used  with  utmost  care,  and  if  any,  Chauvenet's  Criterion,  provides 
the  best  results  without  degradation  in  the  ability  to  represent  the  actual  balance 
behaviour. 


8  CALIB  -  The  Computer  Program 

A  computer  program,  written  in  the  C  and  X /Mob/ programming  languages  has  been 
developed  for  use  in  the  wind  tunnels  at  AMRL.  The  aim  of  the  computer  program  is 
to  allow  effective  and  efficient  analysis  of  various  balance  calibration  matiiematical 
models,  and  to  enable  real-time  conversion  from  balance  voltage  output  to  the 
corresponding  load  experienced  by  the  balance  during  wind  tunnel  tests. 
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8.1  Program  structure 

A  modular  programming  structure  was  used  to  ensure  a  high  level  of  flexibility  in  the 
code.  The  advantages  of  such  a  structure  ensures  minimum  modification  of  the  existing 
data  acquisition  code  if  a  new  balance  calibration  model  is  integrated  into  the  current 
data  acquisition  system.  For  example,  a  third  order  regression  metiiod  for  balance 
coefficient  calculations  can  be  added  easily  to  the  existing  code  as  an  individual 
function,  and  this  would  not  require  any  modification  to  the  existing  code.  A  flow 
chart  for  the  program  CALIB  is  given  in  Section  8.3. 

All  variables  used  in  the  program  are  stored  in  a  data  structure,  hence  the  code  can  be 
easily  integrated  with  the  current  data  acquisition  system. 


8.2  Program  operation 

The  program  requires  a  balance  calibration  data  input  file  which  must  be  in  the  format 

of  {Hi . He,  Ri . Re}  and  a  balance  design  load  range  input  file,  listing  the 

maximum  design  load  range  for  a  particular  balance  imder  investigation  (see 
Appendix  B.l  and  B.2  for  samples  of  these  data  input  files). 

A  data  output  file  shown  in  Appendix  B.3  is  generated  by  the  program.  It  contains 
details  of  the  calibration  model,  optimisation  details,  the  balance  calibration 
coefficients,  standard  errors  and  the  coefficient  of  multiple  correlation. 

The  graphical  user  interface  (GUI)  (see  Appendix  C)  provides  a  user  friendly 
environment  for  the  end  users.  The  GUI  clearly  lists  all  the  available  options  in  a  single 
window,  and  users  may  select  the  calibration  equation  and  appropriate  options  for  tihe 
analysis. 
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8.3  CALIB's  Flow  Chart 


li 
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{CALIB's  flaw  chart  contd) 

This  part  of  the  flaw  chart  represents  the  calculation  of  the  calibration  coefficients  and  reverse 
calibration  for  model:  [H]  =  [C][R-H]  (6, 17, 84, 33, 96  coefficients) 
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{CALlB's  flow  chart  contd) 

This  part  of  the  flow  chart  represents  the  calculation  of  the  calibration  co^cients  and  reverse 
calibration  for  model:  [R]  =  [C][H]  (5, 6, 20, 27, 84, 33, 96  coefficients) 


aaussif ) 

Perform  matrix  inversion  using  Gauss- 
Jordan  elimination  with  full  pivoting 
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{CALIB'sflozv  chart  contd) 


T 
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rev  calibi 


Perform  reverse  calibration 
using  [C]  with  non-optimised 
calibration  data  for 
[H]=[C][R-H] 


rev  calib2i 


Perfoim  reverse  calibration 
using  [C]  with  non-optimised 
calibration  data  for 
[R]=[C][H] 


rev  calibS, 


Perfonn  reverse  calibration 
using  [C]  with  non-optimised 
calibration  data  for 
[H]=[C][R] 


statisticf ) 

Perform  statistical  analysis  on 
the  accuracy  of  the  calibration 


coeff^ ) 

Print  out  all  the  results  to  an 
output  file 


Calibration  data  output 
file 


STOP 
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9  Conclusion 


All  three  calibration  models  have  very  similar  behaviour  in  terms  of  accuracy.  As  the 
order  of  the  calibration  equation  increases,  so  too  does  the  accuracy  of  the  estimated 
forces  and  moments.  A  small,  but  insignificant  improvement  in  accuracy  was  fotmd 
between  the  second  order  (27  coefficients)  and  third  order  (33  coefficients)  calibration 
equations.  The  small  difference  between  these  second  and  third  order  results  may 
primarily  be  due  to  backgrovmd  interference  rather  than  actual  balance  component 
interaction.  A  significant  improvement  in  accuracy  was  found  by  using  the  2"^  order  84 
and  3'^^  order  96  coefficient  balance  calibration  equations. 

Although  a  significant  improvement  in  accuracy  can  be  achieved  by  various  data 
optimisation  techniques,  it  is  recommended  that  data  optimisation  techniques  which 
require  the  elimination  of  calibration  data  should  not  be  used.  This  is  because  those 
data  points  being  eliminated  may  actually  represent  the  true  behaviour  of  the  balance. 
The  results  obtained  from  the  optimised  calibration  matrix  with  the  non-optimised 
calibration  data  set  show  a  reduction  in  accuracy  of  the  estimated  loads.  This  is  due  to 
the  inability  of  the  optimised  calibration  matrix  to  estimate  the  loads  for  those 
calibration  data  points  which  were  eliminated  during  the  optimisation  process.  Of  the 
optimisation  techniques  presented,  Chauvenet's  Criterion  is  the  most  suitable  as  it 
provides  a  reduction  in  standard  error  by  elkninattng  a  minimum  number  of  data 
points. 

To  avoid  the  elimination  of  calibration  data  points  whilst  aiming  to  achieve  a  high  level 
of  accuracy  for  the  calibration  model,  it  is  suggested  that  the  calibration  matrix  should 
only  cover  a  sub-divided  load  range  within  the  balance  design  load  range.  This  would 
lead  to  the  use  of  more  than  one  balance  calibration  matrix  to  cover  the  required  load 
range  of  the  balance. 

Due  to  the  nature  of  the  least  squares  regression  methods,  an  adequate  number  of 
calibration  data  points  should  be  provided  for  the  calculation  of  the  calibration 
coefficients.  It  was  found  that  the  reverse  calibration  method  might  fail  to  converge  if 
an  inadequate  number  of  calibration  data  points  were  provided. 

The  computer  program,  CALIB,  allows  an  efficient  and  effective  way  to  apply  various 
calibration  equations  described  in  this  report  with  different  combinations  of  data 
optimisation  techniques.  This  flexibility  allows  the  user  to  select  the  most  appropriate 
order  of  the  calibration  equation  and  data  optimisation  technique  based  on  the  test 
requirements. 


39 


DSTO-TR-0857 


10  Acknowledgements 

The  authors  would  like  to  thank  the  following  individuals  and  indicate  their  area  of 
contribution. 

Dr  Stephen  Lam  -  software  programming  in  the  C  and  X/Motz/ programming 

languages; 

-  optimisation  technique  for  strain  gauge  balance  calibration. 

Dr  Neil  Matheson  -  providing  comments  and  advice  for  this  investigation. 


11  References 

AIAA,  1995,  Assessment  of  wind  tunnel  data  uncertainty,  American  Institute  of 
Aeronautics  of  Astronautics.  (AIAA  S-071-1995) 

Cook,  T.  A.,  1959,  A  note  on  the  calibration  of  strain  gauge  balances  for  mnd  tunnel  models, 
Royal  Aircraft  Establishment,  Bedford,  England. 

Galway,  R.  D.,  1980,  A  comparison  of  methods  for  calibration  and  use  of  multi-component 
strain  gauge  zvind  tunnel  balances,  NRC  No.  18227,  National  Aeronautical  Establishment, 
National  Research  Coxmcil,  Canada. 

lAI  Engineering  Division,  1998,  SK9514r-08  Balance  Calibration  Matrix,  LAI  Engineering 
Division,  Tel  Aviv,  Israel. 

Lam,  S.  S.  W.,  1989,  A  Fortran  program  for  the  calculation  of  the  calibration  co^cients  of  a 
six-component  strain  gauge  balance,  ARL-TM-410,  AR-005-598,  Department  of  Defence, 
Australia. 

Sprent,  P.,  1969,  Models  in  regression  and  related  topics,  Methuen  &  Co  Ltd,  London. 

Vetterling,  W.  T.,  Teukolsky,  S.  A.,  Press,  W.  H.  and  Flannery,  B.  P.,  1988,  Numerical 
Recipes  in  C,  Cambridge  University  Press. 


40 


DSTO-TR-0857 


Appendix  A  Balance  Calibration  Models 

A  complete  list  of  all  the  balance  calibration  mathematical  models  are  listed  below, 
where  i  =  1 ...  6  (rmless  otherwise  specified). 

A.1  Calibration  Model:  [R]  =  [C][H] 

A.1.1  First  order,  6  component  equation  with  6  coefficients 

R.  =  + c,.  2/?2 + 

A.1.2  First  order,  5  component  equation  with  5  coefficients 

R,  = 


where,  i  =  2...5 

A.1.3  Second  order,  6  component  equation  with  27  coefficients 

/?, = q,H, + 

+  +C,2,H2^  +C,33//3^ 

+ 

+  ^23^2^3  24/f2H4  +C,.  25^2^5  +C,.  26H2/f6 

^1,45^4^5  46^4^6 

+  C,56HsHe 


A.1.4  Second  order,  5  component  equation  with  20  coefficients 

/?,  = 

+ 

Q, 24^2^4  "*■  Q,25^2^5  ^i,26^2^6 

+  C,34//3^4  +C,35i:^3^5  +C,,eH,H, 

+  C, 45^4^5+^45^4//, 

where,  i  =  2...5 


41 


DSTO-TR-0857 


A.1.5  Third  order,  6  component  equation  with  33  coefficients 


+  +C,33//3^  +C,33/// 

^i.23^2^3  ■*■  ^1,24^2^4  Q,25^2^5  26 ^2^6 

+ 

+  +  c, 333^^3^  + 
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A.1.6  Second  order,  6  component  equation  with  84  coefficients 


/?,  =  C,,/f ,  +  +  C,,H,  + 

+ 

■*"  Q,12^1^2  ^1,13^1^3  Q,14'^1^4  ^1,15^1^5 

■*"  ^/,23^ 2^3  ■*"  Q, 24^2^4  ^i.25^2^5  “*■  ^i, 26^2^6 

+  +  C,334^3if  3  +  C, 3,^/3/?, 


•'i,35'‘ 


■*■^1,45^4^5  ■*‘Q,46'^4^6 
+  C,56^5^6 


+q|,3||H,H3|+C,j,3||H,H,|+C,,,.||H,H.|  +  C,j,,||H,H3|+C,,,,||H,H.| 

^/,|23l  1^2^41  ■*‘^ij2S|l^2^5|'*'^J.|2e|  1^2^61 

+  C,j„||H,H4|+C,j„||H2H5|+C,j„||if3H,| 

+  C,,4,||W,ff3|  +  C,|,„|H4H4| 

■*■^<.12611^5^61 

■*■  ^<M|2|^i!^2|'*'^<,||3|^1  l^sl  ■*■^<,1141^1 1^4]  ■*■^<,1151^1  l^s|^*^^<.l|6j^l|^6i 
+  C,_,,H3|H,|  +  C,,3|3|H,|H4KC,.3|3|«3|H5KC,,2|6|H2|«6| 

'*'^<,3|4|^3|^41'*‘^<,3|5|^3|^5|'*'^<J|6|^3|^6| 

+c„|3|H4//,|+c„|^h,|h,| 

■*^^<,3|6|^il^6l 

■K^,j,|2|H,|»2  ■K^,(,|.|ff.|ff4  ■K;,«5|ff,|«5  ■^C,,|„6|H.|^^, 

+  C;i,i,lHi|/r,  +C  ;  loU  ;  lok 


i,|2|3 


•''■.24 


.|2|5l 


+  C,  |3|,|H3lH,  +  C,|3|3|i/3|//3  +  C,|3|g|//3|//g 

+  C,|4|3|//4l//3+C,|,|g|//,|/fg 


+c. 


,|5|6 
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A.1.7  Third  order,  6  component  equation  with  96  coefficients 

~  ^1,2^2  Q,3^3  "*■  ^i,5^5  ^i,6^6 

^i,|l|  1^1 1  ■*■  ^i,\2\  1^2 1  +  ^/,|3|  1^3 1  +  ^i,|4|  1-^4 1  +  ^i,|5|  |  +  ^/,|6|  1^6 1 

+  +  C,33^r3^  +  +  C,esHe" 

+  +  C,3,if  3^ ,  +  +  C,3,^r3^6 

■*■  ^1,34^3^4  ^i,35^3^5  ^,36^3^6 

^/,45^4^5  ■*■  Q,46^4^6 

+  C,,^,,|H,ff,|+C,^,3||//,J^3|+C,^,^|^/,^/,|  +  C,J^|^,H,|+C,,,^ 

*^|J23|  1^2^31  ■*■  ^ij24||^2^4|  ■*'^,J25|  1^2^51  ■*'^/.|26l  1^2^61 

+  CM,4||^?3-ff4|+C,j3„|H3H,|  +  C,,|33||H3H,| 

+  C,j„l|H,A/3i+C,j3,||H,H.| 

+  q,|2|H,|ff2l+C,3|3|H,|ff3|+C,-„ff,|/?.|+C,,,|,H,|fl3|+C,,„fl,|H,| 

'^^i.2|3l^2[^3|'*'^i.2|4]^2|^4|'^^i,2lS|^2|^5!'*'^i,2|6|^2|^6| 

'*'  ^i.3|4l^3 1^4 1  ^j,3|S|^3  i^5 1  '*'  ^i.3l6|  ^3 1^6 1 

■*■  ^i,4|5|'^4|^5|  ■*■  ^(,4|6|-^4|-^6| 

■*"^/.5|6|'^S  1-^6 1 

+  C,J,|2|H,|ff2  +C,„|3|ff,|^f3  +C,,|,|4|«,|H4  +C.^p\H,^  S+ ^ 

+  C,tp\Hz\H,  +  C,(,.|H2|ff4  +  C,h,|H2|H3  + 

+  C,J3|4|»3|»4  +  C,j3|5|H3|Hi  +  C,(3,Jff  jI^  3 

+  C,„<.\»s\H6 

+  C,,,„H,’  +C,„,H,^  +C,^H’  +C„„H,’C,^+H,^ 

,|222|  kl+c, 

,|333|  kl  +  C, 

,|444|  |^4l  +  Qj555||^5^|  +  ^/  ,1666||^6  I 
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A.2  Calibration  Model:  [H]  =  [C1[R] 

A.2.1  First  order,  6  component  equation  with  6  coefficients 

+  Q.2^2  +  +  Q.4^4  +  ^i,5^5  +  ^i,6^6 

A.2.2  First  order,  5  component  equation  with  5  coefficients 

Hj  =  C;  2-^2  3-^3  4-^4  S'^S 


where,  f  =  2... 5 

A.2.3  Second  order,  6  component  equation  with  27  coefficients 


/•  —  ^1,1 -^1  ^1,2^2  Q,3-^3  Q,4^4  ^/,5^5  "*■ 

■^^1,22^2  ■^Q.33'^3  +C,-,44i?4  +  C;  sjiJj  +C;_g6/?5 

^1,12^1^2  Q, 13^1^3  Qm4^1^4  ^1,15^1^5  "*■  ^i,16-^1^6 

Q,23^2^3  ■*"  ^1,24^2^4  ■*"  25^2^5  ■*"  26^2^6 

Q,34'^3^4  ^i,35^3^5  ^/,36^3^6 

Q,45^4^5  ■*■  Q.46^4^6 
■*■  Q,56^5^6 


A.2.4  Second  order,  5  component  equation  with  20  coefficients 

^1  ~  ^/,2^2  ■*■  ^i.3^3  ■*■  ^iA^4  ^1,5^5  ^i.6^6 

■*■^1,22^2  ■*'^j,33^3  ■*■^1,44^4  "^'Q.SS^S  ^i,66^6 

^1,23^2^3  ^1,24^2^4  '*'  ^1,25^2^5  ^i,26^2^6 

^1,34^3^4  Q, 35^3^5  ^i, 36^3^6 

^,45^4^5  ■*■  ^i.46^4^6 
+  ^.56^5^6 


where,  i  =  2...5 
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A.2.5  Third  order,  6  component  equation  with  33  coefficients 


H,  =  C,,R,  +  C,,/?,  +  C,3i?3  +  C,,R,  +  C,5/?5  +  C,,/?, 

■*■^1,22^2  ■*■^1,33^3  +^,.44^4  Q,55^5 

^1,13^1^3  ^1,14^1^4  ■*"  Q,I5^1^5  ^,16^1^6 

^i, 23^2^3  '*’^1,24^2^4  "^^i, 25^2^5  '*'^(,26^2^6 
■*■  ^1,34^3^4  ■*■  ^1,35^3^5  ^,36^3^6 

■*■  ^i,45-^4^5  ■*■  ^i, 46^4^6 
^1,56^5^6 

■^QmII^:  ■^^j,222-^2  +^-,333^3  +  ^(,444'^4  555-^5  Q, 666^6 
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A.2.6  Second  order,  6  component  equation  with  84  coefficients 

~  ^i,2^2  ^1,3^3  Q.5^5  Q,6^6 

1^1 1  ■*■  ^i,|2|  1-^2 1  ^i,\3\  1^3 1  +  ^i,\4\  1-^4 1  +  ^/,|5|  I^S  |  +  ^i,|6|  |^6 1 

+  C,„/?>C,32/?/  +C,33^3^  +C,35/?3^ 

^i,l|l|  -^1 1^1 1  ^i,2|2l  ^2 1^2  K  ^1, 3(31 -^3 1^3 1  "*■  ^iA\4\  ^4 1■^4 1  ^i,5|5|  ^5 1'^S  K  ^i.6|6|  ^6 1-^6 1 

Q,I2^1^2  ^(,13^1^3  ^,14^1^4  ■*"  Q.lS^l^S  '*’  Q, 16^1^6 

^1,23^2^3  Q,24^2^4  25^2^5  ^/, 26^2^6 

■*■  ^1,34^3^4  ^(,35^3^5  Q,36^3'^6 

+  ^,,45^4^5  +  ^,.46^4^6 

+  <^,.56^5^6 

'*’  ^I,|12|  |■^1■^2|  ■*■  ^i,|l3|  |•^1■^3|  ■*■  ^i,|l4|  |■^l•^4l  ^ijl5|  1-^1  ^sl  "*"  ^i,|l6|  l-^l-^el 
^i,|23|  |■^2^3l  ■*■  ^i,|24|  |■^2■^4|  +  ^i,|25|  |■^2■^5l  ■*"  ^i,l26|  |■^2•^6| 

^i,|34|  !•^3^4|  ■*■  ^,-,|35|  l-^S-^sl  ■*■  ^i,|36|  |•^3•^6i 

+  ^,M45ll^4^S  I +  C';.|46|  1^4^61 
■^^i,|S6ll^5^6l 

■*"  ^/,l|2|  -^1 1'^2 1  ^i,l|3|  -^1 1-^3 1  ■*■  ^i,l|4|'^l  1-^4  ]  ^i,l|5|  1-^5 1  ^i.l|6|  -^1 1'^e  1 

^1, 2131-^2 1-^3 1  ■*■  ^i,2|4|-^2|'^4l  ^i,2|5l^2|'^5l  ■*"  ^;,2j6|'^2|'^6| 

■*■  ^/,3|4|'^3|'^4|  +  ^i,3|S|'^3|-^s|  ■•■  ^/,3|6l'^3|-^6l 
■*■  ^(,4|5|^4|^5K  ^i,4|6| ^41-^61 

+  ^;.5!6|'^5K| 

'*'^1, 11121-^1 1■^2  ■*'^i,|l|3 1^1 1-^3  ■*■  ^1, 11141-^1 1-^4  ■*‘^i.|l|5 1^1 1-^5  "*■  ^i,jl|6 1-^1 1-^6 
'*'^i,|2|3l'^2l-^3  ■^^i,|2|4|^2l^4  ^;,|2|5 1-^2 1^5  ^/,|2|6|'^2|^6 

■^^/,|3|4 1^3 1-^4  ■•■^i,|3|5 1-^3 1-^5  ■*"  ^/,|3|6 1-^3 1-^6 
■*"^i,|4|5|-^4|-^5  ■^^i,|4|6 1-^4 1-^6 
■'■^i.lsle  1^5 1-^6 
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A.2.7  Third  order,  6  component  equation  with  96  coefficients 


H i  —  C,  j/?,  +  ^/,3^3  ^iA^4  Q,5^5  ■*"  ^i,6^6 

+c,j,||«J+c,h|j?,|+c,|3||S3|+c,j^|«,|+c,j^|/!,|+ 

+  +C,,33Ji/  +  C,,33fi3=  +C,,„Ji/  +C,,3,S/ 

^1.2121^21^2!  +  ^i,3|3|^3|^3K 

■*■  Q.I2^1^2  ■*■  13^1^3  ■*■  Q, 14^1^4  "*■  ^i, 15^1^5  ^i.I6^I^6 

■^Q,23^2^3  ■*■^1,24^2^4  '^^1,25^2^5  *^^1,26^2^6 
■*■  ^1,34^3^4  ^(.35^3^5  ^(,36^3^6 

^(,45-^4^5  Q,46^4^6 

■*■  ^(.56^5^6 

^i,|i2|  |■^I■^2 1  ■*■  ^;,ji3|  |•^l•^3 1  ■*■  ^;.|14|  1-^1 -^4 1  ■*■  ^/,|15|  1-^1 -^5 1  "*■  ^i,|i6|  I 

■*■  ^i.|23l  |■^2•^3  I  ■•■  ^/J24|  |■^2•^4 1  ^  ^,  J25|  |■^2■^5 1  +  ^,-,|26|  |■^2•^6 1 

'*■  ^i,|34|  |■^3•^4 1  +  ^/j35|  1-^3 -^5 1  ^i,|36|  \^i^6 1 

■*■  ^i,|45|  1^4 -^5 1  ■'■  ^/J46|  1-^4 -^6 1 
■^Q|56||^5^6| 

■*■  ^i.l|2| |■^2 1  ■*■  ^,-,l|3| -^1  I  ■*■  ^i,lj4| -^1 1-^4 1  ^/,l|5| -^1 1-^5 1  ■*"  ^i,l|6l-^l  1-^6 1 

^/,2|3|  -^2 1-^3 1  +  ^/,2l4|  -^2 1^4 1  +  ^/,2|5|  ^2  i'^S  |  ^i.2\6\  ^2 1-^6 1 

^/,3|4|  -^3 1-^4 1  ■•■  ^/,3|5|  -^3  \^5 1  ^i,3|6|  -^3 1-^6 1 

+  -.4151^41^51  +  C',._4|g|i?4|/?6| 


+  ^/.5|6l-^5r6| 

■*'Q|1|2!'^i|^2  +^/.|1|3|-^i|'^3  ■^^;,|i|4|-^i|^4  +^/,|i|5|^i1-^5  ■*■  ^i,|l|6|-^l|-^6 

^/.|2|3  1-^2 1^3  ^/j2|4i-^2!-^4  ^;.|2|5|-^2|^S  ^i,|2|6 1-^2 1-^6 

■'■  ^i.|3|4|^3 1^4  +^,-,|3|5|-^3|^5  ',13161^31^6 

■’■^ij4|5l^4i^5  +Q|4|6|^4|^6 
■*■^1,15161^51^6 

■*■^1,111^1  ■*"^;,222'^2  ■*^^i,333-^3  ■*■^(,444^4  ■*■  Q-,555-^5  ■*"  ^;,666-^6 
■*■^1,11111^1  ■*■  ^i,|222|  -^2  ■*■  ^(,j333|  ^3  ■*"  ^/,|444|  -^4  ■*■  ^i,|555l -^5  ■*"  ^i,|666|  -^6 
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') 


A.3  Calibration  Model:  [H]  =  [C][R-H] 


A.3.1  First  order,  6  component  equation  with  6  coefficients 

~  ^i,2R^2  ^iAR^A  ■*"  ^i.SR^S  ^i,6R^6 

A.3.2  Second  order,  6  component  equation  with  27  coefficients 

^ i  ~  ^i,lR^i  ^i.2R^2  ^i,2R^3  "*■  ^i,4/f^4  ^i,5R^5  ^i,6R^6 

+C,,„H,"  +C,,„H/' 

-  ) 

-  +  C,,^H,H, ) 

-(c,,„h.//3+c,,„h,hJ 

A.3.3  Third  order,  6  component  equation  with  33  coefficients 

~  ^i,lR^l  ■*■  ^i.2R^2  ^i.3R^3  ^iAR^*  ^i,SR^5  ^i.6R^6 

+C,,33H3*  +C, +C,,„/f/  +C,,„H,VC,,^H3’) 
+C,„H3H,  +C,,„H3ff,  +C,^H,H,) 

+C,,333i//  +  C,,,33H3=  +C,,„H/  +C,«.ff/) 
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A.3.4  Second  order,  6  component  equation  with  84  coefficients 

~  ^iAR^4  Q.5/f^5  ■*"  ^i,6R^6 

-  (c,J//,| +C„|H,|+ C,JH,|  +  C„|H.|  +  +  C,,.|H.|) 

-  ) 

-  (c,,,|,|H,|ff,| + q  ,| +c,,3|3|tf  3|h,| + c,,^.|H4|//,| + 

+C,,„H,H3  +C„,//,H  j 

-  (Cu,«jH3  +  +  C,,  +  C,  ) 

-  {c,,^H,H, + qjjHjff  3  +  q,.//,//, ) 

-  (c,.|23|  c,^,^ \h,h,\+c,Jh,hJ;i 

-  (c,M  |«,ff  .1 + C,„  |H,«s  I  +  c,„^ 

-(c,«|H5ff6|) 

-(q.|2|^^.l^f2l+C,3MH.|ff,l+q,|4|».|H,|  +  q|3|ff,|H3|  +  C,,,|^H,|H3|) 
-fc2|3i^^2|H3|  +  C,.,4|»2|»2|  +  C,,,3|H3|i/,| +  q,3|ff3|ff3|) 
-(c,,3M«.l«4|  +  +  C, 

-(Ci,4l,ff4|H5|  +  C,,4,«|i/4|H6i) 

-(c,.3M«s|«.|) 

-(c,«2|h,|H2+c,j,|,|h.|H3 +q,|4|H,|H, +q,|3|H,|//,  +c,j,|,|H,|ff3) 

-(q|3|3|ff2|ff3+C,j3|4|H2|H.  +  C,j3|5|H2|«5+C„3|,|i/3|H.) 

-(C,fl4l^^3|«4  +C,j3|5|H3|ff3  + 

-(c,Ms|ff4|H5+C,(4|4|«4|»4) 

-(C,|3|6l«5|».) 
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A.3.5  Third  order,  6  component  equation  with  96  coefficients 


Hi  —  Ci  inRi  +  Ci  21(^2  ■*" 

-  +  C,„h;  +  C,„H/  +C,«H/  ) 
-{C„,H,H,  +  C,„H,H,  +C,„H,H.+Cu,H,H,+C„,H,H,) 

-  (c„^3|^^^l«3|  +  C,,,.|H3|H4  +  C«|i|«3|H3|  +  C,,,^H3|»6|) 

-  (Cu|4|ff  3|H.I+ C,,3|,H,|//,|  +  C,^^H,\H,\) 

-(c,.=M»3|H6|) 

-(c,j,i,|tf,|H3  +c,„|3|H,|ff3  +c,„|3|h,|//. +c,^,|3|h,|a;, 

-(CiMjiHjHj  +C,«4|HJH4  +C,j3|,|H3|H,  +C,,3|3|fl'3|H,) 

-(c,„4|H3|'f4  +C,,|„|H,|//5  +C,„|J/f3|H.) 
-(c,„,|H,|H,+C„4|3lH,|ff.) 


-c, 


1,111 


H\  1,221^2  ■*■^1,333^3  +  ^;,444^4  ^1,555^5  "^^1,666 


;,|111| 


H, 


+  C: 


,|222| 


H2  ^i,|333| 


+c 


,|444| 


+c.. 


,|555| 


■'■  ^1^666]  H(, 
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Appendix  B  Program  Input  /  Output  Listings 

B.l  Data  Input  File  Sample  -  Subset  of  the  1886  data  points 


^MO^'-ovoociiHoa  VO  LnmiHo»^oc'io'0<N  c-i  w  <x>  •  •  o  oo  o  vo  .oo<N«7VLnwc^<T»mvo!covo^mo>otf> 
<^fNjr>«-ioo^ovo<Nooooooooooor-i.«-irHOOoo  •oc*‘«>ovo\r*vDtn^mi>^mc'a  VO  in  romoa 

(N  •tMCMCNtHHlC'l .  O  -rHtHTHrHf-lT-l^rHrHrHlNOJ  ‘CararaCJCa 

«0 . *000  0  0000000000000  O . O . 

O  OOOOOO  IP  OOOOOOOOOOOO  OOOOO 

<Noop*^r^ooco<r»o\iHo>miHmfnco-ooflDnoocoinvoiHovinu>tHr'incNP-c^oo^f^*fl'^'oVou*pnj'<j'inc-oco 
o*j»Hcao\rHtHinT-<fncor“r»aocoo6r-c-*p*-eo®rHOoocoflOcoco  •co«^oveorH('im(Noorao\ovinvpvoooiHr^o 
OOVOH'O^OOV.  0V(J\THr4r-iiHfH'«-l«-(iHr-lrHW  ‘i-ili-li-IW^f-IOOOflVr-IOO^OXOO^COOVi-POOOVO^OOCvOp© 

(Hi— 1«-|  . . .  »  *0  •  •  •  •  •  ‘IrHiH  ’tHi-I  »H  rH  iH  tH  t-I  W 

*  **0  •  •  ••  *0000000  0  0001  OOOOOO  •  *0 . O* . 

OOOIOOOOOl  I  I  \  I  I  P  I  I  I  I  I  I  I  I  \  I  OOlOOOOOOOIOOOOOOOO 
I  I  I  t  I  I  I  I  I  I  t  p  I  1  I  1  I  i  p  4  I  p  I  I  i 

r-pvovoooM^^fnoo'roo>c^r>eo^covou) 

^^^iovOfHfloinrgr^*-p*-ioc^f>‘r-r'r-o«OiHt-iOLnvni-*oin^.<n^t^incoovir)for^fo*cf^oj«a'ro^mm 
ofnvoo«H<r)^^-loH•-fr-•ooooOf-l<)o<^o^o^coco®cococoo^r•^oo^‘^^^^noln•^oof'lfoOltH 
>  COOOOrH«HTH<HO*-lt-liHc-l»-Pr-lT-aTHi-4«-IOOOOOOOOO»  -fOfOOCOCanjOINfOOaraeNCNOCNfnTn 
g  OOO  OOOOO  •  . . .  *0000  *000  *000000  *000 

..  . . *000000000  0  00000000  0  0  •  •  *o  *  »  *0 . .  ♦  •  • 

gOOOOOOOOP  OOOOOOOOOOOOIOOO 


1 


I  I 


I  I  I 


I  I 


i  i 


^  T]>r-l®VO^^V&®nOMnvOCV<i4'®ONCOTHVO 

r)in®Hr‘0’^vanaofo®r*r*oo»-ioroinoH«'^»oHt*‘^ovvoj'a^f-vovoH^ojo^®iA«vovo 
•.vor*jmE^vooora®,r'jtHr-p«Hrar>^in^n®r^r*“Or^®o>oovonin<r>or-m«4‘ra.cao^vocaojc4^in® 
eoovooovOfHOt-io  *ca<NC'ioifNrMtMnac>atHTHfH«HH»-p«*^iNi-i®THojTHMracn^tn’^r>c^m,inihv0p**p'O 

OOOOO  ‘HHOO  •  •  »  . . ‘Hf-li-IH^HfHHHOOOOOOO  * 

.  *  .  .  -o  *  '  *100000  0000000000000  *0  . *  *  ■  •  *  'O 

>4000001000  P  I  I  I  X  1  I  I  I  I  I  I  P  I  I  I  I  lOIOOOOOOOOOOO  OOOOO 


I  I  I  t  I 


p  I  p 


I 


I  p  I  I  I  I  {  I  p  t  I 


I  I  I  I 


X  mc»<Ncv®®^rHVDtninint-^fno>r-iov®r-4 

^r-lr-49»r^®inOr-ln^Or-*OOCNr'*OVOvn®mtnmt*^OCT>^r''>;!>r>®®in®^<H^^^®t^r-lrH^®r-(® 

omr^fonaovc^ovov®voinvovovovor*‘Ovo®vovo\o®t*“Vot^ovr-r'r-iOf-irHt-H®oo»-pvor'*ca^'S»mr-ro 

THrar‘«cM(N«H^■lr■lr•l(s^anlr^)c^olfSN^a(So^<7«o^o\ooo^o^  'Ot^®®  ®  ®  ®  ®  ®  ®  ®  in  ui  m  ®  tn  m  tn 
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B.2  Balance  Design  Load  Range  Input  File 


#  This  file  contains  the  balance  loading  range  for  the  Aerotech 

#  balance 

#  All  tinit  are  in  N/Niti 

# 

#  Axial  (X) 

1000 

# 

#  Side  (Y) 

2000 

# 

#  Lift  (Z) 

5000 

# 

#  Roll  (L) 

100 

# 

#  Pitch  (M) 

300 

# 

#  Yaw  (N) 

150 
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B.3  CALIB's  Output  File 


»  This  data  output  file  is  generated  by  program:  calib 

«  Balance  Calibration  model:  AMRL's  first  order  calibration  equation 
0  Model  equation: 

«  R  =  C1*HX  t  C2*HY  +  C3*H2  ♦  C4*H1  +  CS*Hm  ♦  C6*Hn 

((No.  of  calibration  coefficients:  6 

#  Mathematical  regression  model  used:  Multivariable  regression 

#  Tolerance  in  the  original  [H]  measurements:  0.000% 

#  Tolerance  in  the  original  (R)  measurements:  0.000% 

#  'Zero  data  point (s)'  filter  NOT  ACTIVATED 

#  Optimisation  for  outlier(s)  elimination:  NOT  ACTIVATED 

If  Chauvenefs  Criterion  for  outlier(s)  elimination:  NOT  ACTIVATED 

#  The  calibration  data  input  file  contains  1886  data  points. 


»  Balance  Calibration  Coefficients  (Transposed) 


^  HX 

COl  9.SS131564e- 
C02  -l.llS86112e- 

C03  -3.55989620e- 

C04  4.67655411e- 

COS  6.46963845e- 
C06  8. 211382106- 


HY 

03  1.540692S3e- 
OS  6.77280604e- 
05  1.388S7130e- 
04  1.03005363e- 
04  3.04742113e- 
04  -2. 831885246- 


HZ 

06  -1.71699722e-0S 
03  7.10106126O-05 
OS  3.24190447e«03 
03  -5.75349267©-04 
06  2.50076814e-04 
04  le0085229Se-05 


HI 

-3.71443729e-06 

-3.18480774e-05 

6.99404024«-0S 

9.30475144e-02 

6.49532542e-04 

2.84751339€-04 


Hm 

9.26901859e-07 

-4.449414890-06 

-1.64403613e-05 

1.038094036-03 

3.97225368e-02 

8.74790649e-04 


# - - - 

#  Normalised  Balance  calibration  Coefficients  (Transposed) 


#  - 
# 

COl 

C02 

C03 

C04 

COS 

C06 


KX 

l.OOOOOOOOe+00 

-1.16705813e-03 

-3.72322841e-03 

4.89111989e-02 

6.76647304e-02 

8.58812993e-02 


HY 

2.27482l57e-04 

l.OOOOOOOOe-t-OO 

2.05021566e-03 

1.52086687e-01 

4.49949S65e-04 

-4.1812S844e-02 


HZ 

-5.29626097e-03 

2.19039806e-02 

1.00000000e4>00 

-1.77472616e-01 

7.71388596e-02 

3.110896S6e-03 


HI 

-3.99197907e-0S 

-3.42277573e-04 

7.51663307e-04 

l.OOOOOOOOe-fOO 

6.98065442e-03 

3.06027884e-03 


Hm 

2.3334407Se-05 
-1.12012355e-04 
-4,13879944e-04 
2.6133628de-02 
l.OOOOOOOOe+00 
2 .202252736-02 


IP - 

#  Sensitivity  Matrix 

IP - 

II  HX 

COl  9.S6131564e-03 
C02  0.000000006400 

C03  0.000000006400 

C04  0.000000006400 

COS  0.000000006400 
C06  0.000000006400 


HY 

0.000000006400 

6,772806046-03 

0.000000006400 

0.000000006400 

0.000000006400 

0.000000006400 


HZ 

0.000000006400 

0.000000006400 

3.241904476-03 

0.000000006400 

0.000000006400 

0.000000006400 


HI 

0.000000006400 

0.000000006400 

0.000000006400 

9.304751446-02 

0.000000006400 

0.000000006400 


Hm 

0.000000006400 

0.000000006400 

0.000000006400 

0.000000006400 

3.972253686-02 

0.000000006400 


Standard  Errors 

HX 

HY 

HZ 

HI 

0.43587058% 

0.16990011% 

0.13213621% 

0.60054345% 

Hm 

0.07482476% 


« - 

#  Coefficient  of  multiple  correlation 

« - 

*  ax 

0.99993341 


HY 

0.99998672 


HZ 

0.99999345 


HI 

0.99992707 


Hm 

0.99999778 


Hn 

•2.379772306-05 

8.734115636-05 

•6.038883266-06 

1.350767316-02 

1.180803396-04 

7.422144826-02 


Hn 

•3.206313486-04 

1.176764386-03 

'8.136304806-05 

-1.819915056-01 

1.590919366-03 

1.000000006400 


Hn 

0.000000006400 

0.000000006400 

0.000000006400 

0.000000006400 

0.000000006400 

7.422144826-02 


Hn 

0.14332177% 


Hn 

0.99999451 
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Appendix  C  Computer  Program's  Graphical  User 

Interface  (GUI) 


Running  CALIB  provides  the  user  with  the  following  window: 


CALIB  -  Strain  Gauge  Balance  Calibration 


Calibration  Model; 


Data  input  File; 


fPlea:e  select  a  .:alibiaticn  modelj 


Select  Model 


(Default  data  input  file;  load.datj 


[Select  data  input  file 


Result  Output  File:  (Default  CALIB  output  file:  coeff.out)  [Select  data  output  file 

r'  Activate 'zero' data  points  filter  'Zero'filter  range  [5S]:  : 

Regression  model  for  Calibration  Coefficients  Calculation: 

fvlultivariable  Regression  '■  Ramaswamy’s  Least  Square  Method 

r  Introduce  tolerance  into  the  R  and  H  measured  values 
Tolerances  in  measured  R  and  H  values; 

d(R) :  T  dCH) ;  i  ; 


Optimisation  Configuration 
No  Optimisation 
C  1  s.E.  Optimisation 
'  "  2  S.E.  Optimisation  , 

3  S.E.  Optimisation 

EXECUTE 


. Activate  Chauvenet's  Criterion 

[C]  Matrix  Computation  Using  Optimised  Data 
' '  .Activate  Reverse  Optimisation 
■:  De- Activate  Reverse  Optimisation 


QUIT 


W  a  i  T  i  n  g  f  <:■  r  c  ■'  m  m  a  n  q 
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By  dirking  ^  separate  window  containing  all  the  available 

calibration  equations  will  appear  as  shown  below. 


CALIB  (Version  1 ,0) 

CALIB  -  Strain  Gauge  Balance  Calibration 


Calibration  Model;  (Please  select  a  calibration  model) 

Data 


Select  Model 


li 


Optimis. 


I  tans 


[H]  -  [C][R-H]  - 


1 .  [H]=[C][R-H];  First  order  calibration  equation  (6  coeffs.) 

2.  [H]=[C)[R-H]:  Second  order  calibration  equation  (27  coeffs.) 

3.  [H]=[C][R-H]:  Second  order  calibration  equation  (S4  coeffs.) 

4.  fHl==[C][R-H]:  Third  order  calibration  equation  (33  coeffs.) 


5.  [H]=iC3[R . Hj:  Third  order  calibration  equation  (96  coeffs.) 


(R]  =  (C][H]- 


6.  [R]=[C][H]:  First  order  calibration  equation  (6  coeffs.) 


Pick  a  calibration  model 


OK  i 


2  S.E.  Optimisation 

3  S.E.  Optimisation 


'  Activate  Reverse  Optimisation 
#!  De-Activate  Reverse  Optimisation 


EXECUTE 


QUIT 


Waiting  for  command 


After  all  the  selections  have  been  made,  the  program  will  begin  the  analysis  by  clicking 
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After  the  program  has  successfully  finished  the  computation,  a  separate  window  will 
appear  on  top  of  the  original  window,  indicating  the  computation  is  completed. 


CALIB  (Version  1.0) 


CALIB  -  Strain  Gauge  Balance  Calibration 

Calibration  Model:  [Hl=iCl[R-H]:  Thii  d  oide:  'lOAfficienTs  definition  (yo  coeffs.)  iselect  Model 


Data  Input  File; 


(Default  data  inputflle;  load.dat)  iselect  data  inputfile 


Result  Output  File:  (Default  CALIE  output  file;  coeff.out)  jselect  data  output  file 
”  Actiuate 'zero' data  points  filter  'Zero' filter  range  [Ss]:  i 

Regression  model  for  Calibration  Coefficients  Calculation: 
gi  Multivariable  Regression  ("  Ramaswamy's  Least  Square  Method 


r  Introduce  tolerance  into  th<  - - 

Tolerances  in  measured  R  and  f ;  |  Computation  Completed!! 
d(R) ;  I  i  I _ _ . . . . . . - . 


Optimisation  Corifiguration 
S>  No  Optimisation 
1  S.E.  Optimisation 
r  2  S.E,  Optimisation 
•  '  3  S.E.  Optimisation 


T!r??,?net's  Criterion 


[C]  Matrix  Computation  Using  Optimised  Data 
C  Activate  Reverse  Optimisation 
m  De- Activate  Reverse  Optimisation 


Results  ouTpLited  That's  ell  folks!! 
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Appendix  D  Graphs 

D.l  Calibration  Model:  [R]=[C][H]  (1886  calibration  data  set) 


First  order  6  coefficient 
equation 

Second  order  27  coefficient 
equation 

Second  order  84  coefficient 
equation 

“’‘“Third  order  33  coefficient 
equation 

Third  order  96  coefficient 
equation 


a.  No  optimisation 


■  First  order  6  coefficient  equation 


■  Second  order  27  coefficient 
equation 

"Second  order  84  coefficient 
equation  (Divergence) 

"Third  order  33  coefficient 
equation 

"Third  order  96  coefficient 
equation  (Divergence) 


b.  Ise  optimisation 

Figure  10.  Balance  Calibration  Model  [R]=[C][H]  with  1886  data  points 
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standard  Error  (%)  Standard  Error  (%) 
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First  order  6  coefficient 
equation 

Second  order  27  coefficient 
equation 

Second  order  84  coefficient 
equation 

Third  order  33  coefficient 
equation 

Third  order  96  coefficient 
equation 


Balance  Component 


c.  2se  optimisation 


First  order  6  coefficient 
equation 

Second  order  27  coefficient 
equation 

Second  order  84  coefficient 
equation 

■*“  Third  order  33  coefficient 
equation 

“""Third  order  96  coefficient 
equation 


Balance  Component 


d.  3se  optimisation 


Figxire  10.  Balance  Calibration  Model  [R]=[C][H]  with  1886  data  points  (cont'd) 
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HX  HY  HZ  HI  Hm  Hn 


First  order  6  coefficient 
equation 

"^Second  order  27  coefficient 
equation 

“*■  Second  order  84  coefficient 
equation 

Third  order  33  coefficient 
equation 

“""Third  order  96  coefficient 
equation 


Balance's  Component 


e.  Chauvenet's  Criterion  optimisation 


First  order  6  coefficient 
equation 

"""Second  order  27  coefficient 
equation 

"""Second  order  84  coefficient 
equation  (Divergence) 

Third  order  33  coefficient 
equation 

■*"  Third  order  96  coefficient 
equation  (Divergence) 


Balance  Component 


/.  Ise  optimisation  with  non-optimised  calibration  data 


Figure  10.  Balance  Calibration  Model  [R]=[C][H]  with  1886  data  points  (cont’d) 
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First  order  6  coefficient 
equation 

Second  order  27  coefficient 
equation 

Second  order  84  coefficient 
equation 

Third  order  33  coefficient 
equation 

Third  order  96  coefficient 
equation 


Balance  Component 


g.  2se  optimisation  with  non-optimised  calibration  data 


HX  HY  HZ  HI  Hm  Hn 


First  order  6  coefficient 
equation 

■**  Second  order  27  coefficient 
equation 

Second  order  84  coefficient 
equation 

^  Third  order  33  coefficient 
equation 

■*“  Third  order  96  coefficient 
equation 


Balance  Component 


h.  3se  optimisation  with  non-optimised  calibration  data 


Figure  10.  Balance  Calibration  Model  [R]=[C][H]  with  1886  data  points  (cant'd) 
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0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0 

HX  HY  HZ  HI  Hm  Hn 
Balance  Component 

i.  Chauvenet's  Criterion  with  non-optimised  calibration  data 
Figure  10.  Balance  Calibration  Model  [R]=[C][H]  with  1886  data  points  (cont'd) 
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standard  Error  (%)  Standard  Error  [%] 


DSTO-TR-0857 


D.2  Calibration  Model:  [H]=[C][R]  (1886  calibration  data  set) 


Balance  Components 


First  order  6  coefficient 
equation 

Second  order  27  coefficient 
equation 

Second  order  84  coefficient 
equation 

Third  order  33  coefficient 
equation 

“•“Third  order  96  coefficient 
equation 


a.  No  optimisation 


Balance  Component 


First  order  6  coefficient 
equation 

“•"Second  order  27  coefficient 
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Figure  11.  Balance  Calibration  Model,  [H]=[C][R]  with  1886  data  points  (cont'd) 
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Figure  11.  Balance  Calibration  Model,  [H]=[C][RJ  with  1886  data  points  (cont'd) 
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Figure  11.  Balance  Calibration  Model,  [H]=[C][R]  with  1886  data  points  (cont'd) 
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Figure  11.  Balance  Calibration  Model,  [H]=[C][R]  with  1886  data  points  (cont'd) 


67 


standard  Error  (%)  Standard  Error  [%] 


DSTO-TR-0857 


D.3  Calibration  Model:  [H]=[C][R-H]  (1886  calibration  data  set) 


Balance  Components 


First  order  6  coefficient 
equation 

“*■  Second  order  27 
coefficient  equation 

■•“Second  order  84 
coefficient  equation 

■*^  Third  order  33  coefficient 
equation 

■*^  Third  order  96  coefficient 
equation 


a.  No  optimisation 


Balance  Component 


First  order  6  coefficient 
equation 

■•■  Second  order  27  coefficient 
equation 

■•“Second  order  84  coefficient 
equation  (Divergence) 

■^  Third  order  33  coefficient 
equation 

■*■  Third  order  96  coefficient 
equation  (Divergence) 


b.  Ise  optimisation 


Figure  12.  Balance  Calibration  Model,  [H]=[C][R-H]  with  1886  data  points 
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Figure  12.  Balance  Calibration  Model,  [H]=[C][R-H]  with  1886  data  points  (cant'd) 
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Figure  12.  Balance  Calibration  Model,  [H]=[C][R-H]  with  1886  data  points  (cont'd) 


70 


standard  Error  (%)  Standard  Error  (%) 


DSTO-TR-0857 


First  order  6  coefficient 
equation 

■^Second  order  27 
coefficient  equation 

“*■  Second  order  84 
coefficient  equation 

Third  order  33  coefficient 
equation 

■*“  Third  order  96  coefficient 
equation 


g.  2se  optimisation  with  non-optimised  calibration  data 


First  order  6  coefficient 
equation 

Second  order  27 
coefficient  equation 

Second  order  84 
coefficient  equation 

Third  order  33  coefficient 
equation 

Third  order  96  coefficient 
equation 


Balance  Component 


h.  3se  optimisation  with  non-optimised  calibration  data 


Figure  12.  Balance  Calibration  Model,  [HJ=[C][R-H]  with  1886  data  points  (cont'd) 
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Figure  12.  Balance  Calibration  Model,  [H]=[C][R-H]  with  1886  data  points  (cont'd) 
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Figure  13.  Balance  Calibration  Model,  [R]=[C][H]  with  329  data  points 
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Figure  13.  Balance  Calibration  Model,  [R]=[C][H]  with  329  data  points  (cpnt'd) 
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Figure  13.  Balance  Calibration  Model,  [R]=[C][H]  with  329  data  points  (cont'd) 
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D.5  Calibration  Model:  [H]=[C][R]  (329  calibration  data  set) 
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Figure  14.  Balance  Calibration  Model:  [H]=[C][R]  with  329  data  points 
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Figure  14.  Balance  Calibration  Model:  [H]=[C][R]  with  329  data  points  (cont'd) 
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Figure  14.  Balance  Calibration  Model:  [H]=[C][R]  with  329  data  points  (cant'd) 
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D.6  Calibration  Model:  [H]=[C][R-H]  (329  calibration  data  set) 
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Figure  15.  Balance  Calibration  Model:  [H]=[C][R-H]  with  329  data  points 
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Figure  15.  Balance  Calibration  Model:  [H]=[C][R-H]  with  329  data  points  (cant'd) 
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Figure  15.  Balance  Calibration  Model:  [H]=[C][R-H]  with  329  data  points  (cont'd) 


81 


standard  Error  [%] 


DSTO-TR-0857 


D.7  5  Component  Calibration  Model:  [R]=[C][H]  (1886  calibration 

data  set) 


HX  HY  HZ  HI  Hm  Hn 


Balance  Component 

Figure  16.  5  component  calibration  model:  [R]=[C][H]  with  1886  data  points 


82 


standard  Error  [%] 


DSTO-TR-0857 


D.8  5  Component  Calibration  Model:  [H]=[C][R]  (1886  calibration 

data  set) 


Balance  Component 

Figure  17.  5  component  calibration  model:  [H]=IC][R]  with  1886  data  points 
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